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RESUMO 

O objetivo do trabalho foi caracterizar e determinar a variabilidade espacial de atributos 
ffsicos de um Argissolo Amarelo, utilizando metodos geoestatfsticos e confec?ao de 
mapas de superffcie pelo metodo da Krigagem, em duas camadas de solo. Foram 
realizadas coletas de solo em intervalos regulares de 2,5 m, numa area de 400 m 2 , sendo 
81 pontos em cada uma das camadas (0 - 0,2 m e 0,2 - 0,4 m), totalizando 162 amostras 
que foram submetidas as determinagoes ffsicas e os dados obtidos analisados pela 
estatfstica descritiva para verificar a normalidade da distribuigao de freqtiencia e 
posteriormente submetidos a geoestatfstica. O modelo que melhor se ajustou a curva do 
semivariograma, na maioria das variaveis, foi o esferico. No entanto, para a variavel 
densidade das partfculas e agregados com diametro 2,00-1,00 mm, na profundidade de 0 
- 0,2 m e grau de flocula?ao, densidade do solo e umidade obtida com tensao 1 ,5 MPa, 
na profundidade de 0,2 - 0,4 m o modelo que melhor se ajustou foi o exponencial. O 
efeito pepita relativo para o grau de floculagao, densidade do solo e densidade das 
partfculas foi de 100%, caracterizando o efeito pepita puro, confirmando a aleatoriedade 
dos dados. Os maiores alcances foram observados nas profundidades de 0,2 - 0,4 m 
para percentagem de silte (51,00 m) e percentagem de areia (35,06 m). Os valores da 
densidade do solo foram considerados relativamente altos sendo de 1,69 e 1,65 Mg.m' 3 
para 0 - 0,2 e 0,2 - 0,4 m, respectivamente. Houve uma maior retengao de umidade nas 
camadas inferiores. 

Palavras-chave: Geoestatfstica; krigagem; semivariograma; dependencia espacial. 




SUMMARY 



The objective of this research was characterize and determine to spatial variability of 
physical attributes of a Alfisols, using methods geostatistical e and making of surface 
maps for the method of Kriging, in two soil depths. They were made soil collections in 
regular intervals of 2.5 m, in an area of 400 m 2 , being 81 points in each one of the layers 
(0 - 0.2 m and 0.2 - 0.4 m), to add 162 samples that were submitted to the physical 
determinations and the obtained data analyzed by the descriptive statistics to verify the 
normality of the frequency distribution and later on submitted to the geostatistical. The 
model that better it was adjusted to the curve of the semivariogram, in most of the 
variables, it was the spherical. However, for the variable density of the particles and 
aggregate with diameter 2.00-1.00 mm, in the depth of 0 - 0.2 m and flocculation 
degree, bulk density and soil-water content with tension 1.5 MPa, in the depth of 0.2 - 
0.4 m the most appropriate model that better adjusted it was the exponential. The nugget 
ratio for the flocculation degree, bulk density and density of the particles was of 100%, 
characterizing “random variance”, confirming the independence of the data. The largest 
ranges were observed in the depths of 0.2 - 0.4 m for lime percentage (51.00 m) and 
percentage of sand (35.06 m). The values of the bulk density were considered relatively 
high being of 1.69 and 1.65 Mg.m' for 0 - 0.2 and 0.2 - 0.4 m, respectively. There was 
a larger humidity retention in the inferior layers. 

Key words: Geostatistic; kriging; semivariogram; space dependence. 
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1. INTRODUgAO 



Diante de uma crescente demanda de alimentos por parte do mercado consumidor 
globalizado e exigencias cada vez maiores por produtos que zelem pela qualidade, 
produtores estao cada vez mais conscientes de seu papel na agricultura. Buscam 
medidas para otimizar a produgao com a diminuigao das perdas, aumentando a 
potencialidade produtiva do solo. 

Com essa finalidade, e imprescindfvel a adogao de medidas que visem o aumento 
da eficacia dos insumos, incluindo sua aplicagao, alem de utilizar a terra de maneira 
racional e sustentavel, conservando e melhorando progressivamente seus recursos 
naturais. O manejo adequado auxilia na manutengao de caracterfsticas e propriedades 
ffsicas e qufmicas do solo, pois, estes, sao fatores que estao ligados diretamente com o 
aumento e a manutengao da produtividade. Afora isso, os mercados mundiais, alem da 
qualidade dos produtos agrfcolas, passaram a exigir controles sobre todo o sistema de 
produgao, e isso inclui a forma com que o solo vem sendo manejado, ou seja, as 
tecnologias empregadas para a viabilidade da produgao. 

Portanto, ha uma crescente preocupagao, por parte da sociedade, com a 
sobrevivencia do homem no planeta, o que tem concretizado conceitos como produgao 
sustentavel na qual procura-se adequar a atividade agrfcola a uma agao que seja 
ambientalmente correta, socialmente justa e economicamente viavel (Hansen, 1996). No 
ano de 1927, quando a populagao mundial girava em torno de dois bilhoes de 
habitantes, ja era consenso entre os cientistas que expandir a area cultivada nao seria 
solugao para alimentar toda a populagao; competiu aos estudiosos, notadamente da 
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ciencia do solo, viabilizar meios para aumentar a capacidade produtiva (Keeney, 2000). 
O mesmo autor destaca a importancia desta ciencia na produ?ao de alimentos e seus 
desafios diante de demandas igualmente prioritarias, mas antagonicas. 

A estatfstica e uma ferramenta bastante utilizada no sentido de verificar as 
similaridades e dissimilaridades dos atributos do solo. Yieira & De Maria (1994), citam 
a necessidade de se fazer uso de grandes parcelas e com muitas repeti?oes, para se ter 
maior sensibilidade a heterogeneidade do solo. No entanto, devido as exigencias da 
estatfstica classica, proposta por R. A. Fisher, dificilmente serem atendidas em condigoes 
de campo, ja que a maioria das propriedades do solo varia no espago e no tempo, devido 
seu carater anisotropico, tornou-se necessario ferramentas estatfsticas adicionais que 
considerem essas variagoes; a geoestatfstica tem se mostrada bastante oportuna nestas 
situagoes. 

Deste modo, a aplicagao das tecnicas de geoestatfstica permite analise de 
dependencia espacial, como tambem, estimativa de dados, para locais nao amostrados, 
atraves de um estimador sem influencia de fatores externos, para isso se faz necessario 
detalhar a area a ser analisada e conhecer o posicionamento geografico de cada ponto 
amostrado ou o primeiro ponto, caso os demais sejam eqiiidistantes, pois cada atributo 
do solo tera suas informagoes associadas a um dado ponto no espago, permitindo a 
interpolagao, bem como a projegao dos resultados com base na estrutura de sua 
variabilidade natural. A utilizagao da geoestatfstica apresenta grande importancia 
economica, pois quando se conhece a dependencia dos dados e possfvel estabelecer, de 
maneira mais adequada, a implantagao da malha amostral. 

Neste sentido, a hipotese de trabalho foi de que pode haver variabilidade espacial 
de atributos ffsicos do solo, originada tanto por fatores de formagao do solo quanto por 
seu contmuo manejo; trabalha-se, tambem, com a hipotese de que ha uma variancia 
finita dentro do dommio espacial estipulado. 

O presente trabalho foi conduzido em um solo classificado como Argissolo 
Amarelo Eutrofico abruptico, textura arenosa/media, A moderado, mesoeutrofica 
(. Alfisol ), com o objetivo principal de avaliar a variabilidade espacial de alguns atributos 
ffsicos do solo, usando metodos geoestatfsticos, a fim de fornecer subsfdios para o 
planejamento de experimentos conduzidos na area, auxiliando na escolha dos pontos 
amostrais, amparado pela teoria das variaveis regionalizadas. 
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2. REVISAO DE LITERATURA 



2.1. O solo e a agricultura 

Segundo a Sindipetro, a agricultura hoje ocupa 25% da superffcie terrestre, sendo 
que destes, 75% esta em solo considerado de baixa fertilidade, metade em solos 
acidentados, com uma produtividade modesta (Educar, 2005). 

O uso cada vez mais intenso das terras agrfcolas, aliado ao desenvolvimento de 
tecnologias modemas de cultivo e as pressoes economicas para reduzir custos de 
produgao, tem ocasionado mudanga na mentalidade quanto as praticas de cultivo do 
solo. Assim, a correla?ao entre parametros flsicos do solo em fun?ao dos sistemas de 
manejo e do desenvolvimento de rafzes na rizosfera adquire importancia cada vez maior 
(Crestana et al., 1994). 

No entanto, a utiliza?ao de equipamentos pesados, a intensa mobiliza?ao do solo, 
oriunda da sua mecanizagao excessiva e, principalmente, o preparo do solo realizado em 
condigao de umidade inadequada vem causando, em muitas regioes, a degradagao da 
estmtura original, ocasionando sua compactagao, i.e., alteragao no arranjamento das 
partfculas, tornando-o mais denso. Essa compacta?ao altera varias de suas propriedades, 
notadamente as relacionadas aos seus aspectos flsicos, como o aumento da resistencia a 
penetragao de rafzes e da densidade do solo, redugao da porosidade total, principalmente 
a macroporosidade e parametros hfdricos, ocasionando, com isso, alteragoes, na maioria 
das vezes indesejaveis, no espago poroso do solo (Silva et al., 1986). 
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De acordo com Marcos (1968), o estudo da estmturagao do solo no sentido 
morfologico e puramente descritivo, enquanto que do ponto de vista ffsico, e funcional. 
Segundo o autor, embora os agregados sejam separados, sua estabilidade avaliada e sua 
formagao extensamente estudada, nao ha conhecimento de alguma lei ou princfpio 
comandando a sua formagao, tamanho e composigao. Tentativa nesta diregao foi 
reportada, entre outros, no trabalho de Harris et al. (1966). 

Alem do intenso uso do maquinario agrfcola e abrasao por implementos de 
cultivo, Reichert et al. (1993) relacionam outras for?as responsaveis pela desagregagao 
do solo, sao elas: (a) impacto de agua (gotas de chuva) e cisalhamento pela agua 
(escoamento), e (b) entrada de agua nos agregados (Kemper, 1965). O metodo de 
peneiramento em agua simula adequadamente condigoes de irriga?ao em sulcos, na qual 
ocorre umedecimento dos agregados e for?as de cisalhamento (Kemper & Chepil, 
1965), mas nao simula a desagregagao por impacto de gotas de chuva (Glanville & 
Smith, 1988). 

Silva & Mielniczuk (1997) ressaltam que a formagao e a estabiliza?ao da estrutura 
do solo ocorrem mediante a atuagao de processos fisicos, qufmicos e biologicos. Esses, 
por sua vez, atuam por mecanismos proprios, nos quais sao envolvidas substancias que 
agem na agregagao e na estabilizagao. Entre essas, as principais sao: argila, sflica 
coloidal, compostos organicos, metais polivalentes, carbonato de calcio, oxidos e 
hidroxidos de ferro e alummio, exsudados organicos e substancias organicas 
provenientes da agao dos microrganismos. Alem das substancias agregantes, existem os 
agentes de agregagao, representados por clima, rafzes, microrganismos e pelo proprio 
tracionamento do solo. 

Segundo Vieira & Vieira (1981), e necessario que se utilize o solo da maneira 
mais racional possfvel, buscando uma metodologia adequada que expresse as 
possibilidades do meio e que represente um aproveitamento equilibrado do ecossistema, 
fazendo-se necessario lembrar que para lograr exito na exploragao do solo e necessario 
deter o conhecimento de suas possibilidades e da relagao que existe entre elas e o meio 
ambiente. 

2.2. Variabilidade espacial dos atributos fisicos do solo 

A preocupagao com a variabilidade espacial remonta ao infcio do seculo passado, 
quando Smith (1910) estudou a disposigao de parcelas no campo em experimentos de 
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rendimento de variedades de milho, numa tentativa de minimizar o efeito de variagoes 
no solo. Montgomery (1913), preocupado com o efeito do nitrogenio no rendimento do 
trigo, fez um experimento de 224 parcelas nas quais mediram o rendimento dos graos. 

Segundo Bacchi & Reichardt (1992), a variabilidade espacial de atributos ffsicos 
dos solos e decorrente da agao diferenciada dos processos de genese e de manejo. 
Quando se trata da ciencia do solo, as dificuldades tornam-se maiores, porque, sendo a 
formagao do solo resultante de uma serie de fatores (fatores de formagao do solo), o 
produto final e bastante heterogeneo em relagao as suas caracterfsticas e propriedades. 
O assunto tem sido abordado em seus diversos aspectos e bastante explorado nos 
ultimos anos, atraves da ferramenta da geoestatfstica. 

Cruz (2005) relata que a heterogeneidade e inerente ao solo, o que tipifica sua 
anisotropia. Em consequencia, o solo apresenta ampla variabilidade espacial de seus 
atributos, tanto no sentido espacial quanto no temporal, a qual se deve ao material de 
origem diverso e as condigoes variadas de sua formagao e pode ser aleatoria ou 
depender da posigao das amostras no espago. Segundo Souza (1992), e importante deter 
o conhecimento da variagao espacial de variaveis ffsicas e qufmicas do solo para 
levantamentos e classificagao de solos, avaliagao da fertilidade, desenvolvimento de 
esquemas de amostragem mais adequados, para o planejamento experimental, bem 
como para as definigoes de praticas apropriadas de manejo e de recuperagao. 

No sentido de praticas apropriadas de manejo, Grego & Vieira (2005) acrescentam 
a importancia da variabilidade espacial na pesquisa agrfcola, que estuda o solo e a sua 
capacidade produtiva, pois a consideram o princfpio basico para o manejo preciso das 
areas agrfcolas, qualquer que seja sua escala. Nesse contexto, Souza (1992) aponta que a 
variabilidade pode ainda ser alterada pelo manejo do solo. 

As propriedades hfdricas tambem sao afetadas pela anisotropia do solo. Segundo 
Warrick & Nielsen (1980) e Webster & Olivier (1992), a condutividade hidraulica e 
uma das propriedades do solo que apresenta maior variabilidade espacial. Por isso 
recomenda-se um grande numero de observagoes para uma analise relativamente 
confiavel (Van Lier & Libardi, 1999). Outros trabalhos abordam a variabilidade 
espacial das propriedades hfdricas do solo, como os de Vieira et al. (1981) e Logsdon & 
Jaynes (1996). 

Timm (2002) explica que a variabilidade espacial do solo pode ocorrer em 
diferentes nfveis, relacionados a diferentes fatores, tais como: variagao do material de 
origem, clima, relevo, organismos e tempo, i.e., relacionado com os processos de 
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formagao do solo e efeitos das praticas de manejo adotadas no uso agrfcolas. 
Ferramentas estatfsticas autocorrelogramas, crosscorrelogramas, semivariogramas, 
analise espectral, krigagem, co-krigagem, modelos autoregressivos, modelos 
autoregressivos integrado media movel (ARIMA), modelos de espago de estados etc., 
tem sido usados para estudar a variabilidade espacial dos atributos do solo, e podem, 
potencialmente, conduzir a praticas de manejo que permitam um melhor entendimento 
do processo de interagao no sistema solo-planta-atmosfera. 

Segundo Chichota et al. (2003), acredita-se que as variagoes aumentem com a 
distancia entre amostras e que uma parcela explicavel pela dependencia espacial esteja 
embutida na variagao geral da propriedade ou caracterfstica. Assim, o valor de 
determinado atributo numa posigao poderia ser predito, em parte, pelo valor do mesmo 
atributo numa posigao proxima. Tecnicas como o autocorrelograma podem ser 
utilizadas para testar esta hipotese. 

2.3. Estatfstica classica, geoestatfstica e a ciencia do solo 

Em 1956, Fisher propos uma teoria de analise estatfstica segundo a qual o 
problema da variagao poderia ser resolvido atraves da casualizagao e repetigao, e tem 
sido referida como “ estatistica classica” (Nielsen et al., 1983). Ela considera o campo 
em termos de conjunto de observagoes, reduzidas a uma media e medidas da sua 
divergencia, sem considerar a sua distribuigao espacial. Para atender tais exigencias, 
procura-se selecionar um local tlpico e representativo, onde as amostras sao retiradas ao 
acaso, evitando-se locais espacialmente variaveis. Na maioria das vezes, entretanto, 
essas exigencias nao sao cumpridas, havendo comprometimento na avaliagao dos 
resultados. Devido as limitagoes da estatlstica classica, de suas exigencias nunca terem 
sido testadas em condigoes de campo e serem dificilmente satisfeitas, ja que a maioria 
das propriedades do solo varia no espago e no tempo, tornou-se necessaria utilizagao de 
ferramentas estatlsticas adicionais que considerem e reflitam essas variagoes. 
Alternativamente, pode-se empregar a teoria das variaveis regionalizadas, proposta por 
Matheron (1963, 1971), para mineragao e posteriormente para ciencia do solo (Vieira et 
al., 1981; McBratney & Webster, 1983). A caracterfstica basica dessa teoria esta no fato 
de se utilizar dessa estrutura da variancia, no espago ou no tempo, para avaliagao do 
domlnio representado por uma amostra. Na utilizagao da geoestatlstica, ha necessidade 
de conhecer a localizagao relativa das amostras no espago ou no tempo, a partir do que 
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se pode calcular a estrutura da variancia, a qual expressa a dependencia entre as 
medidas, atraves de autocorrelogramas e semivariogramas. Assim, pode-se determinar a 
distancia ou o tempo de decorrencia entre as amostras, o qual representa o dommio da 
medida. Quando os dados apresentarem periodicidade ou houver concentragao de 
variancia em determinados periodos, a analise espectral pode ser utilizada para 
determina-los (Silva et al., 1989). 

Ate recentemente, as pesquisas na area agronomica concentravam-se no uso da 
estatfstica classica (analise de variancia, media, coeficiente de variagao, analise de 
regressao etc.), que pressupoe a independencia entre as observagoes, desconsiderando o 
local de amostragem no campo. Neste caso, os experimentos sao conduzidos tendo em 
vista a minimizagao do impacto da variabilidade espacial, isto e, ignorando o fato de 
que as observagoes podem ser dependentes espacialmente. Nielsen & Alemi (1989) 
comentam que o fato das observagoes dentro e entre tratamentos nao serem 
independentes entre si, em muitos casos, toma o desenho experimental inadequado. Tais 
observagoes podem apresentar graus de variabilidade bastante altos, o que dificulta a 
analise e ate pode prejudicar o experimento de tal forma que a parcela seja considerada 
inadequada para a experimentagao feita. Estas variagoes podem, no entanto, estar 
relacionadas com a variabilidade espacial das caracterfsticas ou propriedades do solo, 
tanto vertical como horizontalmente, mesmo em parcelas aparentemente homogeneas 
(Cadima et al., 1980; Libardi et al., 1986; Reichardt et al., 1986; Albuquerque et al., 
1996). 

Para evitar o problema de analisar o grau de dependencia espacial, sao 
empregados, muitas vezes, metodos como a utilizagao de areas consideradas 
homogeneas e a adogao de repetigoes distribuldas no terreno ao acaso. Neste caso, 
sendo verificada a distribuigao normal de freqiiencias, assume-se, sem testes adequados, 
a independencia entre as amostras (Vieira, 2000; Chichota et al., 2003). 

Embora seja mais conhecida em problemas de geologia e mineragao, a teoria das 
variaveis regionalizadas, comumente conhecida como geoestatlstica, tambem pode ser 
aplicada a outras ciencias, como por exemplo, a do solo, promovendo a obtengao da 
estrutura da variancia, porque considera a dependencia entre as medidas, o que nao 
ocorre na estatlstica classica. Assuntos especfficos de real interesse e importancia 
agronomica incluem: umidade, textura, densidade, disponibilidades de agua as plantas, 
condutividade hidraulica, velocidade de infiltragao, elevagao, precipitagao, temperatura 
(Libardi et al., 1986). 
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A analise geoestatfstica e usada para definir a dependencia espacial das 
propriedades do solo tanto para a isotropia quanto anisotropia (Burgess & Webster, 
1980). Propriedade da superficie do solo (Brejda et al., 2000), conteudo de nutrientes 
dos solos (Lee et al., 2001), lixiviagao de nitrato (Ersahis, 2001), distribuigao de 
pesticidas no solo e parametros ecologicos (Rossi et al., 1992) sao exemplos de 
pesquisas analisadas pela geoestatlstica para prognosticar a variagao espacial das 
propriedades do solo. 

O entendimento da variabilidade exige uma maior quantidade de informagoes, as 
quais podem ser obtidas a partir de operagoes de amostragem de solo no campo. Quanto 
maior o numero de amostras, mais alto e o onus financeiro para o agricultor, e, portanto, 
o metodo da geoestatfstica e uma alternativa na otimiza?ao deste processo (Mulla et al., 
1992). 

Conhecido o semivariograma da variavel, e havendo dependencia espacial entre as 
amostras, podem-se interpolar valores em qualquer posigao no campo de estudo, sem 
tendencia e com variancia mmima. Esta interpola?ao chama-se krigagem, em 
homenagem ao engenheiro de minas D.G. Krige (Vieira, 2000). 

Trabalhando com dados de concentragao de ouro, Krige (1951) concluiu que 
somente a informagao dada pela variancia e insuficiente para explicar o fenomeno em 
estudo. Para tal, e necessario levar em consideragao a distancia entre as observagoes. A 
partir daf surgiu o conceito da geoestatfstica, que leva em consideragao a localizagao 
geografica e a dependencia espacial (Ortiz, 2002). 

A geoestatfstica se difundiu bastante na Ciencia do Solo durante a decada de 80, 
quando diversos trabalhos ajudaram a expandir essa ferramenta, como os realizados por 
Reichardt (1984) e Vieira et al. (1983). Outros autores apresentam a geoestatfstica como 
ferramenta para o estudo de propriedades do solo e determinagao do numero de 
amostras a serem retiradas de uma area (Warrick et al., 1986; Trangmar et al., 1985; 
Silva, 1998 e Ortiz, 2002). 

Segundo Gongalves et al. (1999), a teoria das variaveis regionalizadas permite a 
modelagem da dependencia espacial de uma propriedade, podendo ser usada para 
estimar valores medios com mmima variancia por meio da krigagem. A variancia da 
estimativa depende somente do grau de dependencia espacial, ou seja, do 
semivariograma e da configuragao dos pontos vizinhos usados na krigagem. 

Ortiz (2002) considera que o estudo da dependencia espacial de variaveis do solo, 
atraves da teoria das variaveis regionalizadas, permite a interpretagao e a projegao dos 
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resultados com base na estrutura da sua variabilidade natural. O estudo da variabilidade 
de um solo, alem de caracterizar uma regiao, pode indicar o numero e a distribui?ao de 
amostras a serem retiradas possibilitando ainda maior detalhamento da area e dos 
resultados. 

Webster & Olivier (1990) afirmam que muitas propriedades do solo variam 
continuamente no espago e, conseqiientemente, os valores em locais mais proximos 
entre si tendem a ser mais semelhantes, ate um determinado limite, correspondente ao 
dommio destas propriedades, que aqueles tomados a maiores distancias. Caso isto 
ocorra, os dados nao podem ser tratados como independentes e um tratamento 
estatfstico mais adequado e necessario, pois o tratamento estatfstico puro da 
heterogeneidade ignora a existencia da correlagao espacial entre os valores medidos. A 
analise geoestatistica permite descrever esta correlagao, o que torna possfvel, em 
aplicagao pratica, a interpolagao por krigagem (Vieira et al., 1983), alem de permitir a 
identificagao de padroes de amostragem adequados (McBratney & Webstney, 1983; 
Fonteh & Podmore, 1994). 

Por isso, Vieira et al. (1983), McBratney & Webster (1983) e Souza et al. (1997), 
preconizam que o estudo da dependencia espacial das propriedades do solo pode reduzir 
o numero de amostras em relagao ao uso dos procedimentos amostrais definidos na 
estatfstica classica. Com esta preocupagao em mente, a amostragem casual simples 
poderia ser evitada, visto que existe a probabilidade de coletar amostras em locais 
similares, duplicando dessa forma a informagao. 

Ainda neste contexto, Libardi et al. (1986) descrevem que as coletas de amostras 
na estatlstica classica sao feitas ao acaso, ignorando-se a sua posigao geografica. Na 
geoestatlstica, os locais de amostragem sao importantes porque se considera a 
dependencia entre as medidas. Nesse sentido, a geoestatfstica e uma ferramenta que nao 
pode ser ignorada, uma vez que e capaz de propiciar o desenvolvimento de tecnologia 
para a amostragem e descrever a variabilidade das propriedades do meio flsico de um 
sistema, alem de possibilitar a interpretagao dos resultados com base na estrutura da 
variabilidade natural dos atributos avaliados, considerando a dependencia espacial 
dentro do intervalo de amostragem. 

Teixeira et al. (2000), destacam a importancia de se realizar uma coleta bem feita, 
pois quando a coleta de amostras de solo nao e bem-sucedida, todo o trabalho posterior 
fica comprometido, seja por problemas de insuficiencia numerica de amostras, seja pela 
sua baixa qualidade. Para determinar alguns dos parametros edafologicos (umidade, 
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fertilidade, granulometria, entre outros), o uso de amostras deformadas e suficiente, 
mas, em alguns casos, e necessario manter inalteradas as condigoes do solo “in situ”, ou 
seja, em sua estrutura natural, caso em que sao absolutamente necessarias as amostras 
indeformadas. 

Segundo Orlando Filho & Rodella (1983), cerca de 80% a 85% do erro total nos 
resultados usados na recomendagao de fertilizantes e corretivos podem ser atribufdos a 
amostragem no campo e de 15% a 20% podem ser decorrentes do trabalho de 
laboratorio. Daf a necessidade de um bom plano amostral. 

Reichardt et al. (1986) consideraram que para obter dados representativos sobre 
suas caracterfsticas ffsicas, mecanicas e qufmicas, se faz necessario retirar grande 
numero de amostras, o que dificulta a fase operacional da coleta, logo a necessidade de 
uma malha amostral adequada se faz necessaria. 

Segundo Trangmar et al. (1985), a ferramenta geoestatfstica permite descrever, 
com maior precisao, uma variavel qualquer, nao so por diminuir a variagao aleatoria a 
ela atribufda, mas tambem por definir a distancia mmima entre as amostras para que 
possam ser consideradas independentes. Desta forma, a analise da variabilidade e da 
dependencia espacial de uma variavel e importante para subsidiar o planejamento de 
experimentos, bem como para avaliar os efeitos dos tratamentos por reduzir a variagao 
experimental atribmda ao erro aleatorio. 

Conhecendo as coordenadas geograficas do ponto amostrado podem-se analisar os 
dados, possibilitando, segundo Vendrusculo (2001), representar a area com maior 
detalhamento. O nfvel de detalhamento e obtido por meio da distancia entre pontos de 
amostragem e depende tanto da propriedade a ser analisada quanto da escala de 
trabalho, i.e., tamanho da area amostrada (Grego & Vieira, 2005). 

Em experimentos realizados sob condigoes rigidamente controladas, como casa de 
vegetagao e/ou laboratorios, e esperado que o coeficiente de variagao (C.V.) seja baixo, 
geralmente menor que 5%. Mas para experimentos conduzidos sob condigoes onde nao 
se pode ter um controle rfgido o C.V. pode apresentar valores altos. 

O semivariograma e definido como sendo a metade da diferenga do quadrado 
entre valores da amostra numa dada distancia “lag space” (variancia), i.e.: 

/(h) = \ {Var[Z(t) - Z(t + h) ] 2 } (1) 




11 



Segundo Guimaraes (2001), a Var[Z(t) - Z(t + h)] e a variancia dos dados por uma 
distancia h, mas, na expressao acima, esta variancia esta sendo dividida por dois, entao 
se utiliza o prefixo “semi” para distinguir da variancia e dal vem o nome semivariancia 
para y(h) e semivariograma para o grafico de y(h) em fungao de h. 

A caracteriza§ao do semivariograma e dada pelo alcance da dependencia espacial 
(a), que e a distancia limite da dependencia espacial, sendo que as medigoes realizadas a 
distancias maiores que (a), tem distribui§ao aleatoria e, portanto, sao independentes 
entre si. 

O Efeito Pepita (C 0 ) e o valor que revela a descontinuidade do semivariograma 
proximo a origem em distancia menores do que a menor distancia de amostragem 
(Vieira et al., 1983). Quando o Efeito Pepita ( C 0 ) for aproximadamente igual ao Patamar 
(C 0 + C), tem-se o efeito pepita puro, o que demonstra uma ausencia total de 
dependencia espacial (Tragmar et al.,1985). 

Segundo Vieira et al. (1983), diferentes tipos de erros aleatorios ou sistematicos, 
na obtengao das amostras e medidas das variaveis em intervalos inferiores ao de 
amostragem, geram uma descontinuidade na origem, denominada de efeito pepita 
(“ nugget effect”). Um exemplo de semivariograma pode ser visto na Figura 1. 

O Patamar (C 0 + C), corresponde a variancia da amostra (s 2 ) para dados 
estacionarios. O patamar e, portanto, o ponto onde toda semivariancia da amostra e de 
influencia aleatoria. Deste ponto em diante, considera-se que nao existe mais 
dependencia espacial entre as amostras. 




Figura 1. Semivariograma com efeito pepita, Ortiz (2002). 
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A letra minuscula “a” da Figura 1 e denominado de alcance, e a zona de 
influencia da amostra. Ele marca a distancia a partir da qual, para determinada variavel, 
um ponto nao tem influencia sobre outro vizinho, ou seja, inicia-se a aleatoriedade para 
aquele atributo. Esse tipo de situagao do semivariograma refere-se ao crescimento que 
tende ate o patamar, que marca o valor da variancia total da variavel em estudo. 

2.4. Influencia dos atributos do solo na sua produtividade 

Tendo como escopo a manutengao da produtividade do sistema agrfcola, 
condigoes satisfatorias das propriedades fisicas devem ser desenvolvidas, pois somente 
o uso de fertilizantes, melhoramento genetico e medidas para controle de pragas e 
doengas, nao preservarao a produtividade, se ocorrer uma degradagao significativa das 
condifoes fisicas (Oliveira et al., 1983). 

Segundo Carpenedo & Mielniczuk (1990), a desagregagao estrutural e evidenciada 
pelo aumento da densidade do solo, redugao do tamanho dos agregados, da 
macroporosidade, da taxa de infiltragao de agua e do desenvolvimento radicular das 
plantas, varios pesquisadores desenvolveram pesquisas neste sentido, dentre eles 
destacam-se trabalhos desenvolvidos por Cintra (1980); Silva (1980); Machado & Brum 
(1978); D’Agostini (1981) e Dalla Rosa (1981). 

O preparo do solo pela aragao e gradagem e o uso subseqiiente de cultivares sao 
realizados, na maioria das vezes, fora das condi?oes ideais de umidade. 
Conseqiientemente, sao esperadas modifica?6es nas propriedades fisicas, prejudicando a 
capacidade produtiva do solo (Byrnes et al., 1982). As modificagoes que sofrera a 
estrutura serao diferentes de solo para solo e, no mesmo solo, serao diferentes com a 
variagao do equipamento empregado e da intensidade de seu preparo. Os efeitos serao 
na diminuigao da infiltragao das aguas de chuva, aumentando o escorrimento superficial 
e possibilitando maior erosao (Lal, 1979). 

Quando o solo passa de uma condigao natural, mato ou campo nativo, para o 
cultivo anual, ocorrem modificagoes nas propriedades ffsicas, sendo estas mais afetadas 
quanto maior a intensidade de preparo do solo. As principais alteragoes sao 
evidenciadas pela diminuigao de macroporos, tamanho dos agregados, taxa de 
infiltragao de agua e aumento da densidade do solo (Souza & Cogo, 1978; Dalla Rosa, 
1981; Machado et al., 1981; Reinert et al., 1984; Ros et al., 1997). 
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Conforme Kiehl (1979), quanto mais elevada for a densidade do solo maior sera 
seu grau de compacta?ao, menor sera sua estruturagao, menor sua porosidade total e, 
consequentemente, maiores serao as restrigoes para o crescimento e desenvolvimento 
das plantas. Archer e Smith (1972), relatam que a densidade do solo e uma importante 
propriedade para avaliagao das condigoes de campo, alem de ser um elemento 
indispensavel para o calculo da porosidade total do solo, para conversao de expressao de 
conteudo de agua, auxiliar na avaliagao da estrutura, e que mudangas em seu valor 
afetam a agua disponfvel, a capacidade de aeragao do solo, a sua taxa de drenagem e a 
penetragao de rafzes das plantas. Existem, atualmente, tres metodologias para a 
determinagao deste atributo, a saber: metodologia do anel ou cilindro volumetrico, do 
torrao impermeabilizado e da proveta. 

O aumento da densidade implica em efeitos desfavoraveis que causam 
impedimento ffsico no desenvolvimento radicular e restringem o movimento da agua e 
do ar ao longo do perfil, por conseguinte intensificam a erosao dos solos. Todas essas 
restrigoes geram limitagoes no uso e manejo com implicagoes na diminuigao da 
produtividade e impactos negativos no ambiente (Silva et al., 1998). 

2.5. Caracterizagao fisico-hidrica do solo 

Segundo Jorge & Prado (1988), a caracterizagao ffsico-hfdrica do solo e essencial 
para nortear as praticas de manejo do mesmo, destacando-se irrigagao, drenagem, 
cultivos e subsolagem. As duas caracterfsticas fisicas responsaveis pela porosidade que 
atuam de forma mais decisiva na economia da agua no solo sao a textura e a estrutura. O 
conhecimento do movimento de agua dentro do perfil de um solo e de fundamental 
importancia para a solugao de problemas envolvendo irrigagao, drenagem, 
armazenamento e transporte de agua e de nutrientes, bem como o controle da infiltragao 
e escoamento superficial, responsavel pela erosao (Pauletto et al., 1988). 

Kutilek & Nielsen (1994) apontam que amostras perturbadas nao podem ser 
utilizadas para caracterizar as propriedades hidraulicas do solo. Os dados obtidos sao 
inuteis, mesmo aqueles obtidos de amostras de solos arenosos e homogeneos. Na 
pratica, porem, observa-se que toda coleta de solo causa alguma perturbagao em seu 
arranjo poroso. 

Varias tem sido as propostas para medir a retengao da agua no solo, estas 
propriedades sao relativamente dificeis de serem medidas (Dirksen, 1999). Os 
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parametros do modelo muitas vezes tornam-se impossfveis de medissem diretamente, 
devido ao tempo, aos recursos financeiros, a instrumenta?ao, a escala e as restrigoes 
conceituais (Abbaspour et al., 2001). Em procedimentos basicos de laboratorio, a 
retengao da agua e medida em amostras em equilfbrio e a condutividade hidraulica do 
solo nao saturado e medida em amostras onde e estabelecido um fluxo estacionario 
(American Society of Agronomy, 1982). Varios metodos sao sugeridos para inferir 
propriedades hidraulicas do solo de dados obtidos em fluxo nao-estacionarios. 

Naime (2001) relata a importancia dos macroporos nos processos solo-planta- 
rafzes. Tais processos tem movimentado pesquisadores a descrever tamanhos e formas 
(Perret et al., 1999). A definigao de macroporo pode parecer simples. No entanto, sua 
defini?ao se torna nebulosa e ambfgua. Nao ha consenso entre os pesquisadores quanto 
a definigao das faixas de tamanho de poros em fun?ao de sua dimensao. Para classificar 
os poros, a grandeza utilizada e o diametro do cilindro equivalente (ECD), derivado da 
area de seqao transversal assumida como sendo circular. Perret et al. (1999) adotaram a 
definigao proposta por Luxmoore et al. (1990): “O termo macroporo inclui todos os 
poros em um perfil que, geralmente, sao drenados ate a capacidade de campo, com 
diametro equivalente maior ou igual a 1 mm”. 
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3. MATERIAL E METODOS 



3.1. Descri^ao da area e do solo em estudo 

A area experimental fica localizada no Campus do Pici, pertencente a 
Universidade Federal do Ceara e esta a 60 m da estagao meteorologica do Centro de 
Ciencias Agrarias, com uma dimensao de 20 x 20 m. O experimento foi conduzido no 
municfpio de Fortaleza, durante o periodo de agosto a dezembro de 2006. As 
coordenadas geograficas do primeiro ponto sao 3° 44’ 39,7” de latitude S e os 
meridianos de 38°34’51,1” de longitude W.Grw (de Greenwich), numa altitude de 19,6 
m referenciada pelo nfvel do mar. 

De acordo com a classificagao de Wilhelm Koppen (1918), o clima e do tipo Aw’ 
e C 2 W 2 AV; Segundo a classificagao de Thomthwaite (1948). Trata-se de uma regiao 
pertencente ao gmpo de clima tropical chuvoso, com temperatura media do mes mais 
frio maior ou igual a 18°C. Os mdices pluviometricos apresentam uma media anual de 
1.390 mm, sendo a precipitagao do mes mais seco menor que 30 mm, onde a epoca mais 
seca ocorre no inverno e o maximo de chuvas ocorre no outono e a temperatura media 
anual e de 26,5°C. 

O solo foi classificado como um Podzolico Vermelho Amarelo equivalente 
eutrofico abruptico A moderado textura franco argilo-arenosa (Coelho & Mota, 1978), o 
que corresponde a um Argissolo Amarelo Eutrofico abruptico, textura arenosa/media, A 
moderado, mesoeutrofico (PAe) no Sistema Brasileiro de Classifica?ao de Solos - 
SiBCS (Embrapa, 2006), um “Alfisol” pelo Sistema Norte Americano de Classifica?ao 
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de Solos - Soil Taxonomy (Buol et al., 1997), e como um “Lixissolo” pelo Sistema de 
Classificagao da FAO-UNESCO (1997), com declividade media determinada pelo 
clinometro de ABNEY de 3%. A area encontrava-se em pousio coberta apenas por 
plantas ruderais, por exemplo, salsa ( Mimosa pudica L), carrapicho ( Cenchrus 
equinatus), pe-de-galinha ( Eleusine indica ), tiririca ( Cyperus rotundus L.), dentre 
outras. Antes do pousio, a referida area foi cultivada com a cultura do melao ( Cucumis 
melo L.). 

3.2. Metodologia 
3.2.1 Trabalhos no campo 

A area selecionada foi demarcada com piquetes e barbante (Figura 2), de maneira 
que a distancia entre um ponto e outro fosse de 2,5 metros, totalizando 81 pontos em 
uma area com dimensao de 400 m 2 (20 x 20 m) (Figura 3). 

Em cada ponto foram realizadas 
duas coletas de solo, com estrutura 
deformada e indeformada (natural), 
sendo esta ultima realizada com o 
amostrador de Uhland e aneis 
volumetricos de ago inox com 
volume de 50 cm 3 , em duas 
profundidades, sendo a primeira de 0- 
0,2 m (camada que sofre influencia 
do implemento agrfcola) e a segunda 
de 0,2-0, 4 m (camada sem influencia 
do implemento agrfcola), totalizando 162 amostras; as amostras com estrutura natural 
foram utilizadas nas determina§oes de densidade do solo, porosidade total, 
macroporosidade, microporosidade e curva de retengao de umidade. As amostras com 
estrutura deformada foram usadas nas deterininagoes de granulometria, argila dispersa 
em agua (argila natural), grau de floculagao, densidade das partfculas e determinagao de 
umidade no aparelho extrator de Richards; paralelamente, foram realizadas analises de 
agregados, no tocante a sua estabilidade, com torroes previamente separados no campo. 




Figura 2. Demarcacao da area. 
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Figura 3. Croqui da area experimental destacando-se os pontos onde 
foram realizadas as coletas de solo. 



Todas as amostras coletadas em campo foram devidamente identificadas com 
etiquetas e acondicionadas adequadamente para o encaminhamento ao Laboratorio de 
Ffsica do Solo, do Departamento de Ciencias do Solo da Universidade Federal do 
Ceara. 

3.2.2 Trabalhos realizados em laboratorio 

3.2.2. 1 Analises ffsicas 

Toda metodologia utilizada neste trabalho para as determinagoes ffsicas das 
amostras de solo esta baseada na literatura pertinente. 

Densidade do solo 

Para determinagao da densidade do solo, neste trabalho, utilizou-se a metodologia 
do anel volumetrico conforme Blake & Hartge (1986). Apos a aplicagao das tensoes 
usadas para determinagao do espago intersticial do solo, as amostras foram secas em 
estufa a uma temperatura de 105° C ate peso constante, depois foram levadas ao 
dessecador para esfriar e posterior pesagem em balanga com duas casas decimais. 
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A densidade do solo e definida como a relagao entre a massa do solo seco e o 
volume total do solo (volume do espago ocupado pelas partfculas solidas mais o espago 
poroso). 

Foi calculada pela formula: 



ps (Mg/m 3 ) = 



peso do solo seco 105°C(g) 
volume total do solo( cm) 



( 2 ) 



Densidade das particulas 

E definida como a rela?ao entre a massa total das particulas solidas de uma 
determinada amostra de solo e o volume ocupado exclusivamente por essas partfculas. E 
de grande importancia no calculo da porosidade total e da velocidade de sedimentagao. 

A densidade das particulas foi determinada seguindo a metodologia do balao 
volumetrico que consiste na determinagao do volume de alcool etflico gasto para aferir 
um balao de 50 mL com 0,02 kg de terra fina seca ao ar (TFSA) em seu interior (Blake 
& Hartge, 1986). Para tal, foram pesados 0,02 kg de solo, colocado em lata de alummio 
de peso conhecido e levado a estufa a uma temperatura de 105° C ate peso constante. 
Depois, levou-se ao dessecador e em seguida, transferiu-se o conteudo de solo para o 
interior do balao com o auxflio de um funil. Aferiu-se com alcool etflico ate a ausencia 
de bolhas de ar e anotou-se o volume gasto na bureta de 50 mL. 

Foi calculada pela formula: 

3 _ massa do solo seco a 105 o C(x) 

PP ( g m ) - _ ( volurne gasto de alcool(mL)) ' ' 



Porosidade total 

Caracteriza o espago apenas em relagao ao volume total de poros do solo. 
Algumas tecnicas ja foram testadas e usadas na determina§ao da porosidade de varios 
materiais porosos, porem apenas poucas delas foram adequadamente avaliadas e aceitas 
para uso em solo (Danielson & Sutherland, 1986). 

A porosidade total (volume total de poros - VPT) foi calculada usando-se volumes 
da densidade do solo e da densidade das partfculas, atraves da equa§ao proposta por 
Vomocil (1965): 
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a = 1 - ( 



£S 

PP 



) 



Em que: 

a - porosidade total 
ps - densidade do solo 
pp - densidade das particulas 



(4) 



Microporosidade e Macroporosidade 

A microporosidade foi determinada pela mesa de tensao, mediante a aplicagao de 
uma tensao correspondente a 0,006 MPa sobre as amostras saturadas, ate que a agua que 
ocupava os macroporos (<|) > 50 pm) fosse drenada. Em seguida, as amostras foram 
pesadas e levadas a estufa a 105 °C ate peso constante (apos a ultima tensao 0,009 
MPa) e, por diferenga de massa de agua, foi estimado o percentual de microporos. A 
macroporosidade foi calculada indiretamente, pela diferenga entre a porosidade total e 
microporosidade, como se segue: 

Macroporosidade = Porosidade Total(m 3 m~ 3 ) - Microporosidade(m 3 m~ 3 ) (5) 



Umidade obtida atraves da mesa de tensao 

Utilizou-se a mesa de tensao para baixas tensoes. As amostras com estrutura 
natural foram previamente saturadas, por um periodo de 24 horas, para serem levadas a 
mesa de tensao, onde se aplicaram diferentes tensoes, a saber: 0,0009 MPa, 0,006 MPa 
(microporosidade) e 0,009 MPa, para a elabora?ao da curva caracterfstica de retengao de 
umidade do solo. Calculou-se a umidade em base de volume: 



Umidade (0) = 



Volume de agua (m 3 ) 
Volume total da amostra (m 3 ) 



( 6 ) 



Umidade obtida atraves do aparelho extrator de Richards 

Foi determinada com o aparelho extrator de Richards, o qual permitiu a utilizagao 
de tensoes mais elevadas. Para este trabalho utilizaram-se as tensoes de 0,033 MPa, o 
que corresponde a capacidade de campo (CC) e 1,5 MPa o que corresponde ao ponto de 
murcha permanente (PMP). 
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Para esta determinagao foram usadas quatro placas de ceramicas porosas, duas de 
1 Bar e duas de 15 Bar, para as tensoes de 0,033 e 1,5 MPa, respectivamente. As placas 
foram saturadas por um perfodo de 48 horas e em seguida dispuseram-se as amostras de 
solo, com estrutura deformada, em aneis de borracha para saturagao por um periodo de 
24 horas (Figura 4). Levou-se ao extrator para que as tensoes fossem admitidas de 
forma gradativa. 

A exemplo da umidade calculada pelo metodo da mesa de tensao, aqui tambem foi 
calculada com base em volume. 



tal tbl 




Figura 4. Aplicagao das tensoes de 0,033 MPa (a) e 1,5 MPa (b) no aparelho 
extrator de Richards, respectivamente. 



Analise granulometrica e argila dispersa em agua 

E a medida da distribuigao por tamanho das parriculas individuais de uma amostra 
de solo. E de suma importancia, pois indica o risco de erosao, disponibilidade de agua 
para as plantas, o uso economico de adubos, a mecanizagao adequada e qual a cultura a 
ser implantada. 

A analise granulometrica foi determinada pelo metodo da pipeta, descrita por Day 
(1965), que se baseia na velocidade de queda das parriculas e dispersao qufmica com 
solugao normal de hidroxido de sodio (NaOH), com agitagao mecanica de 12.000 rpm 
(“ stirrer ”) e escala textural do Departamento de Agricultura dos Estados Unidos 
(USDA), com dispersao de 0,02 kg de TFSA. 
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As areias foram separadas por peneiramento em tamis de malha de 0,053 mm de 
diametro, segundo a classifica?ao granulometrica do Departamento de Agricultura dos 
Estados Unidos. O silte e argila foram separados por sedimentagao das partfculas com 
diametro igual ou inferior a 0,053 mm, conforme a Lei de Stokes. A argila dispersa em 
agua foi determinada pela mesma metodologia, com excegao do hidroxido de sodio 
(dispersante qufmico). A classificagao textural foi determinada seguindo o sistema 
adotado pela Sociedade Brasileira de Ciencia do Solo (Santos et al., 2005). 



Grau de floculagao 

E a relagao entre a argila naturalmente dispersa e a argila total, obtida apos a 
dispersao. Indica a proporgao da fragao argila que se encontra floculada. 

O grau de flocula?ao foi determinado pela equagao: 



GF = [(argila total (g kg' 1 ) - 



argila dispersa em agua (g kg' 1 ) 
argila total (g kg' 1 ) 



)]• 



•(7) 



Estabilidade dos agregados 

Estabilidade de agregados e a resistencia que os mesmos apresentam referente as 
forgas desintegradoras da agua, ar e manipulagao mecanica do solo. 

A analise de agregados foi determinada por via umida, onde foram realizados 
testes para determinagao da sua estabilidade. Para esta determinagao, adotou-se a 
tecnica descrita por Kemper & Rosenau (1986) e Yoder (1936). Os torroes coletados em 
campo passaram por um processo mecanico de peneiramento, onde passavam por uma 
peneira com 8,00 mm de abertura de malha e coletavam-se os agregados na peneira de 
4,76 mm. Em seguida, os agregados eram pesados em subamostras de 0,025 kg e 
umedecidos com jatos intermitentes de agua sob pressao. Apos um perfodo de 24 h os 
agregados foram levados para o oscilador vertical, tambem conhecido como aparelho de 
Yoder, com freqiiencia de 32 oscilagoes por minuto, composto por dois conjuntos de 
peneiras com malhas 2,0, 1,0, 0,50 e 0,25 mm de abertura nas quais ficam retidos 
agregados estaveis. O teor de agregados retidos em cada peneira foi expresso em cinco 
classes de diametro, a saber: 4,76-2,00 mm, 2,00-1,00 mm, 1,00-0,50 mm, 0,50-0,25 
mm e < 0,25 mm. 
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A estabilidade dos agregados foi determinada pela equagao: 



Agregados estaveis = ( 



a - b 
c - (£b) 



) 



( 8 ) 



Em que: 

a - peso do agregado seco a 105° C 
b - peso das areias secas a 105° C 
c - peso da amostra tomada para analise 

3.2.2.2 Analises estatisticas e geoestatisticas 

Os dados foram submetidos a estatfstica descritiva, onde se determinou os valores 
maximos e mmirnos, media, mediana, desvio padrao, coeficiente de variagao (C.V.), 
numero de variaveis, assimetria e curtose. 

Com o programa estatfstico “ Stat ” (STATISTICA, 1999), aplicou-se o teste de 
normalidade de Kolmogorov-Smirnov (K-S), a 5% de significancia e as variaveis que 
nao apresentaram distribuigao de freqiiencia normal foram realizadas transformagoes 
logaritmicas. 

Para a validagao da normalidade dos dados foram tomadas em consideragao as 
medidas de posigao (media e mediana) e assimetria. Para ratificar a valida?ao foi 
aplicado um teste desenvolvido por Morcillo (2004), no qual e levado em consideragao 
o dexperimentah obtido quando do teste de normalidade e o d cn 'ticc» obtido atraves da 
equagao: 

dcritico = 1,36 (ni + n 2 / m x n 2 ) 1/2 (9) 

Em que, para o numero de amostras do trabalho, considera-se m = n 2 = 100. Se, 
dexperimentai ^ dcntico, os dados nao seguiram distribuigao normal, sendo necessario 
transforma-los e caso d eX perimental < d c rftico os dados apresentaram uma distribuigao 
normal e ja se pode aplicar a geoestatfstica com os dados exatamente com estavam. 

Posteriormente, os dados foram submetidos a geoestatfstica. Para cada variavel, 
foram determinados, modelo sugerido, classe espacial de dependencia, efeito pepita 
(Co), patamar (Co+C), alcance e efeito pepita relativo ( nugget ratio), usando 
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semivariogramas, confeccionados com o auxrlio do software GS+ “Geostatistical for 
Environmental Sciences”, versao 5.0.3 Beta (Robertson, 1998). Com este programa 
foram construrdas a estrutura do semivariograma e todos os seus dados. O efeito pepita 
e o fator observado para o grau de dependencia; quanto mais proximo ele estiver da 
origem, maior sera o grau de dependencia. A classificagao do grau de dependencia 
espacial foi feita com base na razao entre o efeito pepita e o patamar (variancia total), 
sendo considerada forte se a razao for < 0,25, moderada quando > 0,25 e < que 0,75, 
fraca se > 0,75 e < que 1 e independente se a razao for igual a 1 (Guimaraes, 2001). O 
alcance foi utilizado para observar a dependencia espacial das variaveis estudadas. 

O modelo utilizado foi sugerido pelo programa estatfstico baseado na regressao e 
na soma dos quadrados dos resfduos, tragou-se a curva para os quatro modelos (linear, 
exponencial, gaussiano e esferico) e o que apresentou um patamar mais claro e definido 
e um menor erro foi sugerido. 

Em seguida, foi realizada a interpola?ao dos dados para observar a distribuigao da 
variabilidade espacial das variaveis estudadas na area, a este processo da-se o nome de 
krigagem. E a partir dele que se produziu os mapas de superfrcie com o programa 
SURFER 8.0 (Golden Software, 2006). 
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4. RESULTADOS E DISCUSSAO 



4.1. Estatfstica classica (distribuigao de frequencia) 

Para que se possa utilizar as tecnicas, tanto da estatfstica de Fisher como a tecnica 
das variaveis regionalizadas, tambem denominada de geoestatfstica, se faz necessario 
que os dados brutos, depois de tabulados e organizados, apresentem uma distribuigao de 
freqiiencia normal, isto e, que as medidas de posigao media, mediana e moda coincidam 
ou, pelo menos, cheguem a uma proximidade. 

No entanto, nao apenas as medidas de posigao media e mediana foram utilizadas 
para a validagao da hipotese de normalidade dos dados. Para tal validagao, baseou-se 
em Souza et al. (2004a) que apontaram os valores de assimetria como indicativo de 
normalidade. Outro fator observado, quando da aplicagao do teste de normalidade de 
Kolmogorov-Smirnov numa significancia de 0,05, foi o d exper i me ntai 0 ue comparado com 
o d cn tico calculado, reforgou ou rejeitou a hipotese de normalidade. O d cn tico calculado 
foi de 0,192. 

Para a granulometria do solo (Tabela 1 e Tabela 2), as variaveis analisadas 
apresentaram distribuigao normal, com excegao do grau de floculagao na profundidade 
de 0-0,2 m e do silte na profundidade 0,2-0, 4 m, por isso, eles foram transformados 
usando transformagoes logarftmicas antes do calculo da semivariancia. Os fatores 
observados para a determinagao da distribuigao nao normal destes atributos foram o 
valor do d ex perimentai, pai'a o grau de floculagao, que foi de 0,203 e a assimetria e o valor 
do d eX p er i menta i, 6,79 e 0,240, respectivamente, para o silte. 
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O coeficiente de variagao, geralmente expresso em percentagem, da ideia da 
magnitude de variabilidade dos dados, porem nada informa quanto a estrutura de 
dependencia espacial das propriedades do solo, o que so e possfvel atraves de tecnicas 
geoestatfsticas. Para esta pesquisa, a variabilidade dos dados foi classificada de acordo 
com os criterios propostos por Warrick & Nielsen (1980), que consideram os valores do 
coeficiente de variagao entre 12% e 60% como media variabilidade e os valores abaixo 
e acima deste intervalo como baixa e alta variabilidade, respectivamente. Neste sentido, 
as variaveis areia, grau de flocula?ao, densidade do solo, densidade das partfculas, 
porosidade total e microporosidade para ambas profundidades apresentaram um C.V. 
abaixo de 12% (Tabela 1, Tabela 2, Tabela 3 e Tabela 4), sendo considerado, portanto, 
como baixa variabilidade; as demais variaveis apresentaram C.V. entre 12-60%. 

Ao analisarem a variabilidade espacial de atributos flsicos do solo, Nielsen et al. 
(1973) concluiram que, para as variaveis argila, areia e densidade do solo, na camada 
superficial do solo, os dados apresentavam uma distribuigao normal. Segundo os 
mesmos autores, em uma area de 150 ha com 120 pontos de amostragem, foi obtido um 
coeficiente de variagao (C.V.) para argila, silte e areia de 42%, 22% e 17%, 
respectivamente. 

Os dados obtidos para as mesmas variaveis citadas acima obtiveram o mesmo 
comportamento decrescente para a camada superficial. No presente trabalho os C.V. 
para a camada de 0-0,2 m, foram, para argila, silte e areia, de 22,78%, 17,81% e 3,79%, 
respectivamente. Observam-se menores valores quando comparado com os dos autores 
citados, isto pode ser atribufdo a diferen§a de tamanho da area de estudo, quase quatro 
mil vezes maior que a do presente trabalho, ao tipo de solo e as amostras retiradas, pois, 
devido a mecanizagao da area, pode-se ter retido amostras do horizonte superficial 
como de horizontes de transi§ao. 

Observou-se um incremento nos teores de argila de quase 69% da camada 
superficial para a subsuperficial, em valores passou de 126 g kg" 1 para 212 g kg" 1 , nos 
teores de silte de 1,59% (63 g kg _1 para 64 g kg" 1 ) e um decrescimo para variavel areia 
na ordem de 10,7%, passando de 809 g kg' 1 para 722 g kg' 1 . Isso se refletiu na 
classifica§ao textural que passou de franco arenosa para franco argilo-arenosa. Outro 
aspecto interessante sobre as argilas e que no solo ocorra predominancia de argilas do 
tipo 1:1, provavelmente caulinita, pois este argilomineral e o principal constituinte da 
fra§ao argila dos solos dos Tabuleiros (Romero & Lima, 2004). Refor§ando, a rela§ao 
silte/argila, responsavel pelo grau de intemperismo, foi de 0,50 e 0,30, nas 
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profundidades de 0-0,2 e 0,2-0, 4 m, respectivamente, alem da densidade das partfculas 
do solo, muito proxima deste mineral de argila. 

O comportamento das variaveis argila dispersa em agua (natural) e do grau de 
floculagao foram inversamente proporcionais, ou seja, a argila floculada diminuiu com a 
profundidade, logo os valores de argila naturalmente dispersa no solo tenderam a 
aumentar. 

Os histogramas das variaveis granulometricas podem ser observados na Figura 5 e 
Figura 6. 

Tabela 1. Parametros para os atributos fisicos do solo: areia, silte, argila total, argila dispersa em 
agua, grau de floculacao para profundidade de 0-0,2 m. 



ATRIBUTOS 



Parametros 


Areia 


Silte 


Argila Argila 

Dispersa em Agua 

i -i 


Grau de 
Floculagao 












Media 


809,00 


63,00 


126,00 


15,00 


876,00 


Mediana 


815,00 


61,00 


121,00 


12,00 


904,00 


Desvio Padrao 


30,74 


11.35 


28,80 


8,08 


71,32 


C.V. (%) a 


3,79 


17,81 


22,78 


54,28 


8,13 


Mmimos 


720,00 


42,00 


77,00 


01,00 


702,00 


Maximos 


859,00 


98,00 


205,00 


40,00 


992,00 


Curtose 


0,80 


0,69 


0,16 


0,81 


-0,49 


Assimetria 


-1.02 


0,93 


0,66 


1,07 


-0,88 


d b 


0,10 


0,13 


0,07 


0,16 


0,20 


a Coeficiente de 


variagao. 










b d experimental. 












Tabela 2. Parametros para os 


atributos flsicos do solo: areia, 


silte, argila total, argila dispersa em 


agua, grau de floculagao para profundidade de 0,2-0, 4 m. 




- 






ATRIBUTOS 






Parametros 


Areia 


Silte 


Argila 


Argila 


Grau de 










Dispersa em Agua 


Floculagao 




















gKg 






Media 


722,00 


64,00 


212,00 


51,00 


756,00 


Mediana 


737,00 


60,00 


197,00 


48,00 


759,00 


Desvio Padrao 


49.64 


20,94 


43,48 


15,15 


51,22 









27 



Tabela 2 , Cont 









ATRIBUTOS 




Parametros 


Areia 


Silte 


Argila Argila 


Grau de 



Dispersa em Agua Flocula^ao 



gkg ' 1 



C.V. (%) a 


6,87 


32,50 


20,42 


29,38 


6,76 


Mmimos 


573,00 


46,00 


137,00 


13,00 


622,00 


Maximos 


801,00 


233,00 


348,00 


98,00 


639,00 


Curtose 


1,05 


54,70 


0,76 


1,35 


1,80 


Assimetria 


-1.03 


6,79 


0,93 


0,80 


0,36 


d b 


0,12 


0,24 


0,16 


0,11 


0,05 



a Coeficiente de varia^ao. 
b d experimental. 



No tocante a densidade do solo, densidade das partfculas, porosidade total, 
macroporosidade, microporosidade e umidade obtida mediante o aparelho extrator de 
Richards, nas tensoes de 0,033 e 1,5 MPa, todas apresentaram distribuigao de freqiiencia 
normal. Em todos os casos as medidas de posigao foram iguais ou muito proximas 
(Tabelas 3 e 4). 

Foi observado neste trabalho que, a variavel densidade do solo apresentou uma 
distribuigao de freqiiencia normal e um baixo coeficiente de variagao, o mesmo foi 
reportado por Kilig et al. (2004), trabalhando em dois Neossolos Fluvicos na Turquia. 
Cruz (2005), em trabalho conduzido num Argissolo cultivado com cana-de-agucar, 
encontrou o mesmo comportamento para referida variavel. 

Segundo Kiehl (1979), para as condigoes dos solos das regioes tropicais, a 
densidade do solo deve estar entre 1,1 a 1,6 Mg m" 3 . Os dados mostram que, esta 
variavel esta no limite superior, sendo de 1,69 Mg m" 3 para profundidade de 0-0,2 m e 
1,65 Mg m" 3 para profundidade de 0,2-0, 4 m, observa-se ainda, que houve um pequeno 
decrescimo com o incremento da profundidade. Isto pode ser um indicativo do infcio da 
formagao de uma camada endurecida na superficie, provavelmente do intenso 
maquinario agrfcola que trafega naquela area quando do infcio ou final de experimentos 
que la sao conduzidos. 

O mesmo aconteceu com a densidade das partfculas, tendo uma diminuigao de 
2,62 Mg m" 3 para 2,59 Mg m" 3 . Na realidade os valores de densidade das particulas eram 
esperados, pois estao relacionados com a mineralogia do solo e este atributo e bastante 
estatico. Em solos tropicais ha uma predominancia de quartzo e feldspatos e a 
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Figura 5. Histogramas relacionados as propriedades areia, silte, argila total, argila dispersa em agua e 
grau de floculagao na profundidade de 0-0,2 m. 
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Tabela 3. Parametros para os atributos fisicos do solo: densidade do solo, densidade das particulas, 
porosidade total, macroporosidade, microporosidade e umidade 0,033 MPa e 1,5 MPa 
para profundidade de 0-0,2 m. 



ATRIBUTOS 



Densidade Porosidade Umidade 



Parametros 


Solo 
Mg 


Particulas 


Total 


Macro 

3 -3 

- m m 


Micro 


0,033 MPa 

m'^ 


1,5 MPa 













Media 


1,69 


2,62 


0,35 


0,15 


0,19 


0,11 


0,06 


Mediana 


1,70 


2,61 


0,35 


0,15 


0,19 


0,10 


0,06 


Desvio Padrao 


0,06 


0,04 


0,02 


0,03 


0,02 


0,02 


0,01 


C.V. (%) a 


3,68 


1,58 


7,11 


24,85 


11,53 


20,20 


22,75 


Mmimos 


1,44 


2,49 


0,30 


0,08 


0,13 


0,06 


0,03 


Maximos 


1,80 


2,69 


0,44 


0,26 


0,24 


0,19 


0,11 


Curtose 


3,30 


0,27 


2,03 


0,66 


-0,20 


1,20 


1,06 


Assimetria 


-1,40 


-0,52 


1,02 


0,77 


-0,23 


0,44 


1,02 


d b 


0,10 


0,15 


0,07 


0,08 


0,06 


0,07 


0,11 



a Coeficiente de varia^ao. 
b d experimental. 



Tabela 4. Parametros para os atributos fisicos do solo: densidade do solo, densidade das particulas, 
porosidade total, macroporosidade, microporosidade e umidade 0,033 MPa e 1,5 MPa 
para profundidade de 0,2-0, 4 m. 



ATRIBUTOS 



Parametros 


De 

Solo 
Mg 


nsidade 

Particulas 

m _ 3 


Total 


- Porosidade 

Macro Micro 

3 -3 

in' m 


Umidade- 

0,033 MPa 1,5 MPa 

m 3 m ' 3 
















Media 


1,65 


2,59 


0,36 


0,13 


0,23 


00,15 




0,09 


Mediana 


1,65 


2,58 


0,35 


0,12 


0,22 


00,14 




0,09 


Desvio Padrao 


0,06 


0,04 


0,02 


0,03 


0,02 


00,03 




0,02 


C.V. (%) a 


4,01 


1,89 


7,24 


28,98 


11,14 


22,44 




22,36 


Mmimos 


1,45 


2,50 


0,30 


0,06 


0,17 


00,10 




0,05 


Maximos 


1,81 


2,69 


0,45 


0,23 


0,29 


00,23 




0,15 


Curtose 


0,56 


-0,86 


1,58 


-0,37 


-0,39 


-0,4 




-0,13 


Assimetria 


-0,42 


0,06 


0,71 


0,36 


0,33 


00,70 




0,59 


d b 


0,04 


0,11 


0,08 


0,08 


0,07 


0,13 




0,09 



a Coeficiente de varia^ao. 
b d experimental. 
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densidade desses minerais gira em torno de 2,5-2,65 Mg m' 3 (Kiehl, 1979). Em 
decorrencia deste carater estatico, este atributo apresentou o menor C.V., sendo de 
1,58% para camada de 0-0,2 m e 1,89% para a camada de 0,2-0, 4 m. O mesmo 
resultado foi encontrado por Cruz (2005), trabalhando num solo de mesma ordem. 
Afora isso, observou-se que a densidade das partrculas apresentou, devido este carater, o 
menor coeficiente de assimetria de todas as variaveis analisadas, sendo de 0,06 para 
camada de 0,2-0, 4 m (Figura 9). 

A porosidade total se manteve praticamente a mesma nas duas profundidades. No 
entanto, os poros com maior diametro ((|» 50 pm) apresentaram uma diminui?ao em seu 
volume, passaram de 0,15 m 3 m' 3 para 0,13 m 3 m' 3 . Os microporos, ao contrario, 
aumentaram seu volume com o incremento da profundidade, passaram de 0,19 m 3 m' 3 
para 0,23 m 3 m" 3 . Este comportamento pode estar relacionado com o aumento no teor de 
argila, alem de compressoes mecanicas que a area vem recebendo, diminuindo assim 
sua porosidade de aeragao e aumentado, evidentemente, a microporosidade, tanto para 
uma camada como para a outra. 

Com a diminui?ao da macroporosidade e com o aumento da microporosidade na 
camada subsuperficial, observou-se um maior acumulo de agua na camada de 0, 2-0,4 m 
de profundidade, que passou de 0,11 m 3 m" 3 para 0,15 m 3 m" 3 , na tensao da capacidade 
de campo (0,033 MPa) e 0,06 para 0,09 m 3 m' 3 , na tensao do ponto de murcha 
permanente (1,5 MPa). 

Os histogramas para estas variaveis podem ser apreciados na Figura 7, Figura 8, 
Figura 9 e Figura 10. 

Para estabilidade dos agregados, verificou-se que os C.V. foram, em geral, os mais 
elevados se comparado com as demais variaveis estudadas (Tabela 5 e Tabela 6). No 
entanto, sua varia?ao foi classificada como media, tendo seu valor maximo de 52,50%, 
correspondente aos agregados com diametro 0,5-0,25 mm, na profundidade de 0-0,2 m. 
Valores semelhantes foram encontrados por Souza et al. (2004b), que encontraram C.V. 
medios para as mesmas variaveis, concordando com estes resultados estao os 
encontrados por Sa et al. (2000). 
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Figura 7. Histogramas relacionados as propriedades densidade do solo, densidade das particulas, 
porosidade total e macroporosidade na profundidade de 0-0,2 m. 
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Figura 8. Histogramas relacionados as propriedades microporosidade, umidade 0,033 MPa e 
umidade 1,5 MPa na profundidade de 0-0,2 m. 
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Figura 9. Histogramas relacionados as propriedades densidade do solo, densidade das particulas, 
porosidade total e macroporosidade na profundidade de 0,2-0, 4 m. 
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Figura 10. Histogramas relacionados as propriedades microporosidade, umidade 0,033 MPa e 
umidade 1,5 MPa na profundidade de 0,2-0, 4 m. 
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Cruz (2005), em sua estatfstica descritiva, observou que os agregados de mesmo 
diametro tambem apresentaram o C.V. mais elevado, em torno de 106,1%, sendo 
classificado como alta variabilidade. 

Os agregados que apresentaram maior estabilidade aos testes submetidos (via 
umida) foram aqueles com diametro maior, ou seja, 4,76-2,00 mm, que foram de 549 g 
kg agregados" 1 e 566 g kg agregados" 1 para as camadas de 0-0,2 e 0,2-0, 4 m, 
respectivamente, o que era esperado em se tratando de um solo com textura franca 
arenosa. O mesmo resultado foi obtido por Flores (2005), trabalhando em um Argissolo 
Vermelho no municfpio de Pelotas (RS), onde os agregados que apresentaram maior 
estabilidade foram aqueles com diametro 4,76-2,00 mm e menor de 0,25 mm, e com 
menor estabilidade aqueles com diametro 0,5-0,25 mm sendo de 47 g kg agregados" 1 e 
49 g kg agregados" 1 para profundidade 0-0,2 m e 0, 2-0,4 m, respectivamente. 

Nota-se que a estabilidade dos agregados apresentou os maiores valores na 
camada subsuperficial, se comparada com a camada superficial. Possivelmente, o 
manejo do solo possa estar contribuindo para tal. 

Os resultados referentes a analise descritiva encontram-se nas Tabela 5 e Tabela 6, 
onde se observa um ajuste normal a distribuigao de freqiiencia para todas as variaveis. 

Os histogramas para as variaveis de analise de agregados podem ser vistos na 
Figura 11 e Figural2. 
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Tabela 5. Parametros para os atributos fisicos do solo: estabilidade de agregados 4,76-2,00mm, 
2,00-1,00 mm, 1,00-0,5 mm, 0,5-0,25 mm e < 0,25 mm para profundidade de 0-0,2 m. 



ATRIBUTOS 



Agregados (mm) 

Parametros 4,76-2,00 2,00-1,00 1,00-0,5 0,5-0,25 < 0,25 

gkg agregados 1 



Media 


549,00 


52,00 


49,00 


47,00 


300,00 


Mediana 


575,00 


46,00 


45,00 


42,00 


279,00 


Desvio Padrao 


154,24 


23,16 


25,04 


24,98 


97,02 


C.V. (%) a 


28,05 


44,25 


50,09 


52,50 


32,33 


Mmimos 


128,00 


5,68 


4,00 


5,00 


104,00 


Maximos 


880,00 


131,00 


123,00 


118,00 


615,00 


Curtose 


0,15 


1,99 


0,59 


0,36 


1,00 


Assimetria 


-0,61 


1,13 


0,87 


0,86 


0,86 


d b 


0,11 


0,11 


0,10 


0,11 


0,11 



a Coeficiente de variagao. 
b d experimental. 



Tabela 6. Parametros para os atributos fisicos do solo: estabilidade de agregados 4,76-2,00mm, 
2,00-1,00 mm, 1,00-0,5 mm, 0,5-0,25 mm e < 0,25 mm para profundidade de 0,2-0, 4 m. 



ATRIBUTOS 



Parametros 


4,76-2,00 


2,00-1,00 


— Agregados (mm) 

1,00-0,5 

- - g Kg agregados 


0,5-0,25 

■l 


<0,25 


Media 


566,00 


68,00 


73,00 


49,00 


242,00 


Mediana 


584,00 


68,00 


68,00 


42,00 


233,00 


Desvio Padrao 


161,38 


27,04 


37,70 


25,84 


85,44 


C.V. (%) a 


28,50 


39,72 


51,07 


51,97 


35,28 


Mmimos 


243,00 


16,00 


11,00 


11,00 


102,00 


Maximos 


849,00 


132,00 


181,00 


124,00 


434,00 


Curtose 


-1,08 


-0,67 


-0,27 


0,13 


-0,96 


Assimetria 


-0,22 


0,33 


0,58 


0,81 


0,26 


d b 


0,09 


0,10 


0,10 


0,12 


0,07 


a Coeficiente de variagao. 
b d experimental. 
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Figura 11. Histogramas relacionados com a estabilidade dos agregados 4,76-2,00 mm, 2,00-1,00 mm, 
1,00-0,50 mm, 0,50-0,25 mm e < 0,25 mm na profundidade de 0-0,2 m. 
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Figura 12. Histogramas relacionados com a estabilidade dos agregados 4,76-2,00 mm, 2,00-1,00 mm, 
1,00-0,50 mm, 0,50-0,25 mm e < 0,25 mm na profundidade de 0,2-0, 4 m. 
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4.2. Geoestatistica (calculo da semivariancia) 

As variaveis do solo mostraram cliferengas na dependencia espacial com 
determinadas analises de semivariancia. Os modelos esferico e exponencial foram os 
que melhor se ajustaram a curva para definir as variaveis, sendo o primeiro mais 
utilizado (Tabela 7). Esta selegao esta de acordo com Carvalho et al. (2003), que 
encontraram estes mesmos modelos para variabilidade espacial de atributos ffsicos do 
solo. 

Para as variaveis areia, silte, argila total, argila dispersa em agua, porosidade total, 
macroporosidade, microporosidade, umidade 0,033 MPa, agregados com diametro 4,76- 



Tabela 7. Parametros dos semivariogramas isotropicos das propriedades fisicas analisadas. 



Atributos 
do Solo 


Classe c 


Modelo 


Efeito Pepita 

c„ 


Patamar 

c„ + c 


Alcance 
A„ (m) 


ER d 

(%) 


Areia (g kg' 1 ) 


0-0,2 m 


F 


Esferico 


1,0000 


16230 


27,42 


0,06 


0,2-0, 4 m 


M 


Esferico 


1033,0 


40040 


35,06 


25,79 


Silte (g kg' 1 ) 


0-0,2 m 


F 


Esferico 


23,700 


132,00 


5,25 


17,95 


0,2-0, 4 m 


F 


Esferico 


0,0184 


0,0832 


51,00 


22,11 


Argila total (g kg' 1 ) 


0-0,2 m 


F 


Esferico 


80,000 


1174,0 


21,64 


6,81 


0,2-0, 4 m 


F 


Esferico 


754,00 


3093,0 


33,98 


24,37 


Argila disp. Agua (g kg' 1 ) 


0-0,2 m 


F 


Esferico 


6,4000 


65,380 


3,15 


9,78 


0,2-0, 4 m 


M 


Esferico 


85,600 


251,30 


11,77 


34,06 


Grau de floculagao (g kg' 1 ) 


0-0,2 m 


I 


EPP 


5073,4 


5073,4 


- 


100,00 


0,2-0, 4 m 


F 


Exponencial 


552,00 


2795,0 


2,75 


19,74 


Densidade do solo (Mg m' 3 ) 


0-0,2 m 


I 


EPP 


0,0039 


0,0039 


- 


100,00 


0,2-0, 4 m 


F 


Exponencial 


0,0009 


0,0044 


0,87 


20,45 


Densidade partfculas (Mg m 
0-0,2 m 


- 3 ) 

F 


Exponencial 


0,0004 


0,0018 


1,63 


22,22 


0,2-0, 4 m 


I 


EPP 


0,0024 


0,0024 


- 


100,00 


Porosidade total (m 3 m' 3 ) 


0-0,2 m 


F 


Esferico 


0,0001 


0,0006 


3,78 


16,66 


0,2-0, 4 m 


F 


Esferico 


0,0001 


0,0007 


2,98 


14,28 


Macroporosidade (m 3 m' 3 ) 


0-0,2 m 


M 


Esferico 


0,0007 


0,0018 


16,69 


38,88 


0,2-0, 4 m 


M 


Esferico 


0,0006 


0,0017 


16,69 


35,29 


Microporosidade (m 3 m' 3 ) 


0-0,2 m 


M 


Esferico 


0,0002 


0,0006 


17,54 


33,33 


0,2-0, 4 m 


F 


Esferico 


0,0002 


0,0010 


27,83 


20,00 


Umidade 0,033 MPa (m 3 m' 
0-0,2 m 


3 ) 

M 


Esferico 


0,0002 


0,0007 


20,21 


28,57 


0,2-0, 4 m 


M 


Esferico 


0,0007 


0,0013 


13,44 


53,84 
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Tabela 7, Cont. 



Propriedade 
do Solo 


Classe 


c Modelo 


Efeito Pepita 

C 0 


Patamar 

c 0 + c 


Alcance 
A 0 (m) 


ER d 

(%) 


Umidade 1,5 MPa (m 3 


m’ 3 ) 












0-0,2 m 


F 


Esferico 


0,0000 


0,0004 


27,94 


0,00 


0,2-0, 4 m 


F 


Exponencial 


0,0001 


0,0008 


14,03 


12,50 


Agregados 4,76-2,00 mm (g kg' 1 ) 












0-0,2 m 


M 


Esferico 


13280 


33300 


30,50 


39,87 


0,2-0, 4 m 


F 


Esferico 


42400 


26090 


2,98 


16,25 


Agregados 2,00-1,00 mm (g kg' 1 ) 












0-0,2 m 


F 


Exponencial 


127 


555,5 


2,12 


22,86 


0,2-0, 4 m 


F 


Esferico 


94 


735,1 


2,98 


12,78 


Agregados 1,00-0,5 mm (g kg' 1 ) 












0-0,2 m 


M 


Esferico 


337 


674,1 


11,94 


49,99 


0,2-0, 4 m 


F 


Esferico 


226 


1430 


2,98 


15,80 


Agregados 0,5-0,25 mm (g kg' 1 ) 












0-0,2 m 


M 


Esferico 


398 


823,9 


29,60 


48,30 


0,2-0, 4 m 


F 


Esferico 


89 


674,9 


2,98 


13,18 


Agregados < 0,25 mm 


(g kg' 1 ) 












0-0,2 m 


F 


Esferico 


4080 


18580 


45,91 


21,95 


0,2-0, 4 m 


F 


Esferico 


1380 


7245 


2,98 


19,04 



c Classe Espacial: F, forte dependencia espacial; M, moderada dependencia espacial; I, ausencia de 
dependencia espacial (efeito pepita puro). 
d Efeito Pepita Relativo “Nugget ratio” - [Co/C 0 +C] x 100. 

Para o grau de floculacao (profundidade 0-0,2 m) e para o silte (profundidade 0,2-0, 4 m), foram 
realizadas transforma^oes logarftmicas; (log (g Kg' 1 )). 
e EPP = Efeito Pepita Puro. 



2.00 mm, 1,00-0,5 mm, 0,5-0,25 mm e menor que 0,25 mm, em ambas profundidades, 
umidade obtida numa tensao de 1,5 MPa, na profundidade de 0-0,2 m e agregados 
2,00-1,00 mm, na profundidade de 0, 2-0,4 m, o modelo que melhor se ajustou foi o 
esferico. Ja para as variaveis densidade das partfculas e agregados com diametro 2,00- 

1.00 mm, na profundidade de 0-0,2 m e grau de floculagao, densidade do solo, umidade 
obtida na tensao de 1,5 MPa, na profundidade de 0,2-0, 4 m, apresentaram um melhor 
ajuste com o modelo exponencial. Isto esta em acordo com Kili? et al. (2004) que, 
estudando a variabilidade de atributos ffsicos, encontraram para o silte e para a 
densidade do solo o mesmo modelo encontrado no presente trabalho. Neste sentido, 
Souza et al. (2004b), encontraram para as variaveis porosidade total, macro e 
microporosidade um melhor ajuste com o modelo esferico, resultados semelhantes 
foram encontrados por Sousa (1999) para as variaveis areia, silte e argila, para ambas 
profundidades. Alem das umidades obtidas a 0,033 e 1,5 MPa, sendo esta ultima para 
profundidade de 0-0,2 m. 

As variaveis grau de flocula?ao e densidade do solo, para 0-0,2 m de profundidade 
e densidade das partfculas na profundidade de 0,2-0, 4 m, apresentaram uma relagao 
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efeito pepita-patamar de 100%, significa dizer que, os valores do efeito pepita e do 
patamar coincidiram, caracterizando assim o fenomeno denominado de efeito pepita 
puro (EPP), onde o alcance e nulo e nao ha dependencia espacial entre as amostras, 
podendo-se inferir que a variabilidade e maxima para estas variaveis e que os pontos sao 
independentes entre si, atendendo o requisito basico para a aplicagao da estatrstica 
classica que, nesse caso, pode ser utilizada sem restri?6es. 

Para Cambardella et al. (1994), o efeito pepita reflete a variabilidade nao 
explicada ou varia?oes nao detectadas pela distancia insuficiente da malha de 

amostragem. Esse parametro pode ser expresso como porcentagem do patamar, com o 
objetivo de facilitar a comparagao do grau de dependencia espacial das variaveis em 
estudo (Trangmar et al.,1985). 

A umidade na tensao de 1,5 MPa, para primeira profundidade apresentou um 
efeito pepita de 0%, indicando que o erro experimental e praticamente nulo e que nao 
existe variagao significante a distancias menores que a amostrada, ou seja, nao ha 
variabilidade espacial para distancias menores que 2,5 m. Ainda sobre a variavel 
umidade 0,033 e 1,5 MPa, observou-se que o efeito pepita foi mais elevado na 
profundidade de 0,2-0, 4 m, indicando que, para esta profundidade ocorre uma 

variabilidade espacial para distancias menores. 

Todas as variaveis apresentaram um grau de dependencia espacial forte (relagao 
efeito pepita/patamar < 25%), i.e., o efeito pepita esta proximo a origem, moderado 
(relagao efeito pepita/patamar > 25% e < 75%) e independente (relagao efeito 
pepita/patamar igual a 100%). Nas quais, para silte, argila total, porosidade total, 
umidade 1,5 MPa, agregados 2,00-1,00 mm e agregados menores que 0,25 mm 

apresentaram dependencia espacial forte para ambas profundidades e areia, argila 

dispersa em agua, densidade das partrculas, apresentaram a mesma dependencia para a 
profundidade de 0-0,2 m, ainda grau de floculagao, densidade do solo, microporosidade, 
agregados 4,76-2,0 mm, agregados 1,00-0,5 mm e agregados de 0,5-0,25 mm para a 
profundidade de 0,2-0, 4 m. 

Dependencia igual foi verificada por Souza et al. (2004b), quando estudavam 
granulometria do solo. Cruz (2005), obteve o mesmo grau de dependencia para as 
variaveis silte e argila total, numa profundidade de 0-0,2 m. Segundo a mesma autora, a 
macroporosidade e a microporosidade, para camada superficial, apresentaram uma 
moderada dependencia, isto foi constatado neste trabalho. Geralmente, uma forte 
dependencia espacial dos atributos do solo pode estar relacionada aos fatores intrmsecos 
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e a distancia entre os pontos; possivelmente a forte dependencia espacial detectada para 
as variaveis citadas acima seja decorrente de qualquer um dos fatores de formagao do 
solo, tais como o relevo, clima, material de origem, atividade biologica e tempo. 

Em se tratando de um Argissolo e sabendo-se que e um solo autoctone, as 
variaveis granulometricas tem ligagoes pedogeneticas, isto justifica o fato destas 
variaveis, apresentarem uma forte dependencia espacial, salvo o teor de areia na 
profundidade 0,2-0, 4 m que, mesmo assim, ficou muito proximo de forte (25,79%). 

Apresentaram grau de dependencia moderado macroporosidade, umidade na 
tensao de 0,033 MPa, para as duas profundidades estudadas, e microporosidade, 
agregados 4,76-2,00 mm, agregados de 1,00-0,5 e agregados de 0,5-0,25 mm para a 
profundidade de 0-0,2 m, por fim, areia e argila dispersa em agua, na profundidade de 
0, 2-0,4 m. As demais variaveis, com exce?ao do grau de floculagao e densidade do solo 
na profundidade de 0-0,2 m e densidade das partfculas na profundidade subseqiiente, 
apresentaram um grau de dependencia espacial forte. 

Os semivariogramas para as variaveis com moderada dependencia espacial 
ajustaram-se melhor ao modelo esferico. Trabalhando em um Argissolo Vermelho 
Amarelo, Carvalho et al. (2003) obtiveram semivariogramas que apresentaram o melhor 
ajuste com o mesmo modelo para as variaveis macroporosidade e microporosidade, 
concordando com trabalhos conduzidos anteriormente por Salviano et al. (1998) e 
Souza et al. (1997). 

Neste trabalho, foram obtidos para a relagao efeito pepita/patamar, aqui 
denominada de efeito pepita relativo ou “nugget ratio”, os valores extremos desde 0%, 
para a variavel umidade na tensao de 1,5 MPa na profundidade de 0-0,2 m e de 100% 
(EPP), para as variaveis grau de floculagao e densidade do solo, na profundidade de 0- 
0,2 m e densidade das partfculas na profundidade de 0,2-0, 4 m. 

A maior parte das variaveis que apresentaram um alcance maior que a area de 
estudo, i.e., maior que 20 metros, se encontra na camada superficial, isto, 
provavelmente estar relacionado com o maquinario agrfcola, pois pode ocorrer a 
diminuigao da espessura do horizonte superficial e algumas amostras da camada de 0- 
0,2 m pode ser parte de horizonte de transigao ou ate mesmo um horizonte B. As 
variaveis onde se verificou tal ocorrencia foram: areia (0-0,2 e 0,2-0, 4 m), silte (0, 2-0,4 
m), argila total (0-0,2 e 0,2-0, 4 m), microporosidade (0, 2-0,4 m), umidade na tensao de 
0,033 MPa (0-0,2 m), umidade na tensao de 1,5 MPa (0-0,2 m), agregados 4,76-2,00 
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mm (0-0,2 m), agregados 0,5-0,25 mm (0-0,2 m) e agregados menores que 0,25 mm (0- 
0,2 m). 

Os valores do alcance foram maximos para o silte na camada inferior, com o valor 
de 51,00 m, seguido dos agregados menores que 0,25 mm, na camada superficial com 
45,91 m e para areia na camada subsuperficial com 35,06 m. 

Observou-se ainda que, quanto mais estreita a relagao do efeito pepita com o 
patamar, menor tende a ser o alcance e maior tende a ser a variabilidade espacial. 
Assim, com exce?ao das variaveis microporosidade, agregados com diametro 4,76-2,00 
mm, 1,00-0,50 mm e 0,50-0,25 mm na camada de 0-0,2 m, areia e argila dispersa em 
agua na camada de 0,2-0, 4 m e macroporosidade e umidade na tensao 0,033 MPa, para 
ambas camadas, as demais variaveis apresentaram uma estreita rela?ao, destacando-se o 
teor de areia, argila total e argila dispersa em agua, na camada superficial, sendo a 
relagao de 0,0006, 0,06 e 0,09, respectivamente. 

Como citado anteriormente, para a aplicagao de geoestatfstica se faz necessario 
que se tenha o 2° momento da estatfstica que e a variancia, em outras palavras tem que 
se ter um patamar claro e bem definido. Nessa linha de raciocmio, percebe-se a 
importancia do patamar, tambem denominado de soleira ou “sill”. E a partir dele que se 
determina o limite de dependencia e independencia dos dados, ou seja, o alcance. Neste 
sentido, todos os semivariogramas gerados para as variaveis em estudo foram 
isotropicos e apresentaram patamar. O semivariograma isotropico leva em considera?ao 
uma volta de 360° em torno do ponto, ou seja, qualquer amostra que se retire do solo 
dentro de uma zona circular, para os determinados alcances, apresentaram valores 
correlacionados espacialmente. 

Tanto a densidade do solo como a densidade das partfculas tiveram os menores 
alcances de todas as variaveis, sendo 0,87 m para 0,2-0, 4 m de profundidade e 1,63 m 
para 0-0,2 m de profundidade, respectivamente. Isto se deve ao fato do efeito pepita e o 
patamar apresentarem uma estreita rela?ao de 0,0009 e 0,0044 para a densidade do solo 
e 0,0004 e 0,0018 para a densidade das partfculas, como pode ser observado no 
semivariograma. Comportamento semelhante ocorreu para a variavel porosidade total, 
sendo as diferengas efeito pepita/patamar de 0,0005 e 0,0006, e alcance de 3,78 e 2,98 
m, respectivamente para a profundidade 0-0,2 e 0,2-0, 4 m. 

Na Figura 13, Figura 14, Figura 15, Figura 16, Figura 17 e Figura 18, sao 
apresentadas as estruturas dos semivariogramas de todas as variaveis estudadas. 
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Areia (g kg -1 ) prof 0-0,2 m 



Silte (g kg _1 ) prof 0-0,2 m 



Modelo Esferico 

(C 0 = 1; Co+C = 1623; A 0 = 27,42; 
r 2 = 0,986; SQR = 24744) 



0,00 



5,50 



11,00 

DistSincia (m) 



16.50 



22,0C 




Argila Total (g kg -1 ) prof 0-0,2 m 



Argila Dispersa em Agua (g kg 1 ) prof 0-0,2 m 



1202 

S02 



601 

301 



□ 




Modelo Esferico 

(C 0 = 80 ; Co+C = 1174; A 0 = 21,64; 
r 2 = 0,979; SQR = 16968) 



o J ■ ■ ' ' ' ' ' ' ' ' ' * — 

0,00 5,50 11,00 16.50 22,00 

Dastincia (m) 






£ 

c 

y? 



7B.8 



55.9 

35.9 
20,0 




Modelo Esferico 

(C 0 = 6,40; Co+C = 65,38; A 0 = 3,15; 
r 2 = 0,005; SQR = 1359) 



0 



0 



0 , 0 ' 

0,00 5,50 11,00 16,50 22,00 

Ditfancia 



Grau de Floculagao (g kg " ) prof 0-0,2 m 



-3 X 



Densidade do Solo (Mg m ) prof 0-0,2 m 



IE-004 



1434 



Efeito Pepita Puro 
(C 0 = 5073,4171; Co+C = 5073,4171) 



2E-QQ4 



IE-004 



Efeito Pepita Puro 
(C 0 = 0,0039; Co+C = 0,0039) 



q ! 

0,00 



5,50 



11,00 

Distincia 



15,50 



22.00 



E+QDO 1 

0,00 



5,50 11,00 15,50 

Distancia (m) 



22,00 



Figura 13. Semivariogramas isotropicos correspondentes as variaveis areia, silte, argila total, 
argila dispersa em agua, grau de flocula^ao e densidade do solo na profundidade de 0- 
0,2 m. 
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Areia (g kg' 1 ) prof 0,2-0, 4 m 



Silte (g kg' 1 ) prof 0,2-0, 4 m 



3678' 



0.00 



Modelo Esferico 

(C 0 = 1033; Co+C = 4004; A 0 = 35,06; 
r 2 = 0,947; SQR = 216677) 



5.50 11.00 

Distancia (m) 



16.50 



22,00 



)D0' 

0,00 



Modelo Esferico 

(C 0 = 0,0184; Co+C = 0,0832; A 0 = 51,00; 
r 2 = 0,767; SQR = 3,25) 



5,50 11,00 16,50 

DistSncia (m) 



22,00 



Argila Total (g kg' 1 ) prof 0,2-0, 4 m 



Argila Dispersa em Agua (g kg' 1 ) prof 0,2-0, 4 m 



2378 f 



21 58 f 



14391 



719 f 



D* — 

0,00 






0 



Modelo Esferico 

(C 0 = 754 ; Co+C = 3093; A 0 = 33,98; 
r 2 = 0,947; SQR = 141916) 



5,50 



11,00 

Distancia (m) 



16,50 



22.00 



5 137 

E 



V 



0,00 



Modelo Esferico 

(C 0 = 85,6; Co+C = 251,3; A 0 = 11,77; 
r 2 = 0,888; SQR = 1337) 



5,50 11,00 13,50 

Distincia (m) 



Grau de Floculagao (g kg' 1 ) prof 0,2-0, 4 m 



Densidade do Solo (Mg m' 3 ) prof 0,2-0, 4 m 



Modelo Exponencial 
(C 0 = 552,0; Co+C = 2795,0; A 0 = 2,75; 
r 2 = 0,508; SQR = 399505) 



Ql 

0,00 



5,50 



11,00 

DistSncia (m) 



16,50 



22.00 



SD.597E-QQ4 



37.948E-004 



5 25,299 E-004 

I 



12.649E-004 



QO.ODOE-hOQO' 

0,00 






Modelo Exponencial 
(C 0 = 0,0009; Co+C = 0,0044; A 0 = 0,87; 
r 2 = 0,016; SQR = 8,829E-07) 



5,50 



11,00 16,50 

Distanoia (m) 



22,00 



Figura 14. Semivariogramas isotropicos correspondentes as variaveis areia, silte, argila total, 
argila dispersa em agua, grau de floculagao e densidade do solo na profundidade de 
0,2-0, 4 m. 
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.3 

Densidade das Particulas (Mg m ) prof 0-0,2 m 




3 -3 

Porosidade Total (m m ) prof 0-0,2 m 



65,577 E -005 
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Modelo Esferico 




(C„ = 0,0001; Co+C = 0,0006; A 0 = 3,78; 


1 5,669 E -005 
OO.DOOE-hOQO 


r 2 = 0,031; SQR = 2,629E-08) 



0.00 5.50 11,00 15,50 22,00 



Distancia 



3 -3 

Macroporosidade (m m ) prof 0-0,2 m 



3 -3 

Microporosidade (m m ) prof 0-0,2 m 



13,93eE-Q04 



14 . 239 E -004 



94 . 930 E -005 



47 . 465 E -005 



OD.DOOE-hOQO 1 

0.00 



Modelo Esferico 

(C 0 = 0,0007; Co+C = 0,0018; A 0 = 16,69; 
r 2 = 0,972; SQR = 2,006E-08) 



5.50 



11,00 15,50 

Distancia (m) 



22,00 



59,370 E -005 



44 . 9 D 2 E -005 



? 29,935 E- 005 ] 
E 



14 . 967 E -005 



OD.DDOE-hOQQ' 
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Modelo Esferico 

(C 0 = 0,0002; Co+C = 0,0006; A 0 = 17,54; 
r 2 = 0,992; SQR = 9,498E-10) 



5.50 



11,00 15,50 

Distancia <m) 



22,00 



Umidade 0,033 MPa (m 3 m" 3 ) prof 0-0,2 m 



Umidade 1,5 MPa (m 3 m" 3 ) prof 0-0,2 m 



57 . 000 E- 005 ' 



5 D. 250 E -005 



33 . 500 E -005 



16 . 750 E -005 



OD.ODOE-hOQO' 

0.00 




(C 0 = 0,0002; Co+C = 0,0007; A 0 = 20,21; 
r 2 = 0,978; SQR = 3,488E-09) 



5.50 



11,00 15,50 

Distancia <m) 



22,00 




11,00 16 f 50 

DistSncia (m) 



Figura 15. Semivariogramas isotropicos correspondentes as variaveis densidade das particulas, 
porosidade total, macroporosidade, microporosidade, umidade 0,033 MPa e umidade 1,5 
MPa na profundidade de 0-0,2 m. 
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Densidade das Particulas (Mg m' 3 ) prof 0,2-0, 4 m 



Porosidade Total ( m 3 m' 3 ) prof 0,2-0, 4 m 
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70,223 E-005 



OO.ODOE-hOQO' 

0,00 



Efeito Pepita Puro 
(C 0 = 0,0024; Co+C = 0,0024) 
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Macroporosidade ( m 3 m' 3 ) prof 0,2-0, 4 m 



Microporosidade ( m 3 m' 3 ) prof 0,2-0, 4 m 
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0,00 



Modelo Esferico 
d (Co = 0,0002; Co+C = 0,0010; A 0 = 27,83; 
r 2 = 0,957; SQR = 2,131E-08) 
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11,00 15,50 

Distancia (m) 



22,00 



Umidade 0,033 MPa ( m 3 m' 3 ) prof 0,2-0, 4 m 



Umidade 1,5 MPa ( m 3 m' 3 ) prof 0,2-0, 4 m 
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E 
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Modelo Exponencial 



M (C 0 = 0,0001; Co+C = 0,0008; A 0 = 14,03; 
r 2 = 0,970; SQR = 5,256E-09) 



5,50 



11,00 15,50 
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22,00 



Figura 16. Semivariogramas isotropicos correspondentes as variaveis densidade das particulas, 
porosidade total, macroporosidade, microporosidade, umidade 0,033 MPa e umidade 1,5 
MPa na profundidade de 0,2-0, 4 m. 
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Estab. Agregados 4,76-2,00 mm (g kg" ) prof 0-0,2 m 



Estab. Agregados 2,00-1,00 mm (g kg" ) prof 0-0,2 m 
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Modelo Esferico 

(C 0 = 13280; Co+C = 33300; A 0 = 30,50; 
r 2 = 0,877; SQR = 3,012E+07) 



5,50 



11,00 

□istancia (m) 



1 S,50 



22,00 



606 



455 



5 303 



/ 
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Modelo Exponencial 
(C 0 = 127; Co+C = 555,5; A 0 = 2,12; 
r 2 = 0,357; SQR = 19697) 
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Estab. Agregados 1,00-0,50 mm (g kg" 1 ) prof 0-0,2 m 



Estab. Agregados 0,50-0,25 mm (g kg" 1 ) prof 0-0,2 m 




Modelo Esferico 

(C 0 = 337; Co+C = 674,1; A 0 = 11,94; 
r 2 = 0,768; SQR = 12245) 
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Modelo Esferico 

(C 0 = 398; Co+C = 823,9; A 0 = 29,60; 
r 2 = 0,804; SQR = 24470) 
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Modelo Esferico 

(C 0 = 4080; Co+C = 18580; A 0 = 45,91; 
r 2 = 0,901; SQR = 6,778E+06) 
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Figura 17. Semivariogramas isotropicos correspondentes a estabilidade dos agregados com 
diametro 4,76-2,00 mm, , 2,00-1,00 mm, 1,00-0,50 mm, 0,50-0,25 mni e < 0,25 inm na 
profundidade de 0-0,2 m. 
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Modelo Esferico 

(C 0 = 4240; Co+C = 26090; A 0 = 2,98; 
r 2 = 0,069; SQR = 5,468E+07) 
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Modelo Esferico 

(C 0 = 94; Co+C = 735,1; A 0 = 2,98; 
r 2 = 0,008; SQR = 43729) 
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Estab. Agregados 0,50-0,25 mm (g kg' 1 ) prof 0,2-0, 4 m 
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Modelo Esferico 

(C 0 = 226; Co+C = 1430; A 0 = 2,98; 
r 2 = 0,11; SQR = 180963) 
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Modelo Esferico 

(C 0 = 89; Co+C = 674,9; A 0 = 2,98; 
r 2 = 0,164; SQR = 54066) 
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Figura 18. Semivariogramas isotropicos correspondentes a estabilidade dos agregados com 
diametro 4,76-2,00 mm, 2,00-1,00 mm, 1,00-0,50 mm, 0,50-0,25 mm e < 0,25 mm na 
profundidade de 0,2-0, 4 m. 




51 



4.3. Geoestatistica (mapas de superficie) 

Nesta fase da pesquisa, o semivariograma teve um papel primordial, pois alem de 
mostrar a dependencia espacial, serviu como ferramenta basica para a interpola?ao dos 
dados e confecgao dos mapas de superficie que foram modelados a partir do programa 
computacional Surfer 8.0. 

Foi a partir dos mapas que se pode observar o comportamento das variaveis dentro 
da area de estudo. Os mapas de distribuigao espacial das partfculas solidas do solo estao 
apresentados na Figura 19 e Figura 20. 

Nota-se que os mapas das variaveis areia (19a), silte (19b) e argila total (19c) na 
profundidade de 0-0,2 m, apresentaram valores entre 720 e 850 g kg' 1 , 40 e 100 g kg' 1 e 
70 e 210 g kg' 1 , respectivamente (Figura 19). Na profundidade de 0,2-0, 4 m, para as 
mesmas variaveis (20a, 20b e 20c), os valores variaram entre 560 e 820 g kg" 1 , 40 e 230 
g kg" 1 e 130 e 350 g kg' 1 , respectivamente. A variavel areia (0, 2-0,4 m) apresentou a 
maior amplitude, sendo de 260 g kg" 1 . 

A variavel areia teve uma diminuigao com o incremento da profundidade tendo 
seu pico passando de 860 para 800 g kg' 1 , no entanto, obteve os maiores valores em 
ambas profundidades quando comprado com as outras duas variaveis granulometricas 
que seguiram tendencia contraria, aumentando de 80 para 200 g kg" 1 no caso do silte e 
de 200 para 300 g kg' 1 , no caso da argila. Percebe-se que na porgao Sudeste e Nordeste, 
coordenadas 4mx2me 19 m x 2 m, respectivamente, os teores de areia foram mais 
elevados na profundidade de 0-0,2 m e na por?ao Noroeste para profundidade 0,2-0, 4 
m. As variaveis areia e argila total, para as duas profundidades, se comportaram de 
forma inversa, onde nos picos de argila pode-se observar as maiores depressoes nos 
teores de areia. Um fato e que onde os teores de argila foram mais elevados ocorreu 
maior resistencia a penetra?ao do anel volumetrico, devido a diminui?ao do espago 
poroso do solo. 

Para profundidade de 0,2-0, 4 m, observa-se que o silte (20b), apresentou uma 
distribuigao uniforme, exceto nas coordenadas 2 m x 4 m, bem onde se encontrou os 
menores valores de areia e argila total. 

As variaveis argila dispersa em agua (Figura 21a e 22a) e grau de floculagao 
(Figura 21b e 22b) apresentaram uma relagao espacial inversa. Na profundidade de 0- 
0,2 m, para argila naturalmente dispersa e grau de floculagao os valores foram entre 0 e 
40 g kg" 1 e 700 e 1000 g kg" 1 , respectivamente. Na profundidade de 0,2-0, 4 m, os 
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valores variaram entre 15 e 100 g kg' 1 , 620 e 940 g kg' 1 , respectivamente. Os picos nos 
teores de argila naturalmente dispersa encontram-se distribuldos por toda area. No 
entanto, seus maiores valores observados foram onde a argila nao se encontrava 
floculada. Isto pode ser verificado na porgao central tendendo para o Norte, na 
profundidade de 0-0,2 m (Figura 21ab) e na porgao Norte, na profundidade de 0,2-0, 4 m 
(Figura 22 ab). 

A densidade das partlculas, na camada de 0-0,2 m (Figura 21c), foi a variavel que 
mais uniformemente se distribui na area, tendo seus valores variando entre 2,49 Mg m' 3 
e 2,70 Mg m' 3 , para 0-0,2 m de profundidade e 2,50 Mg m' 3 e 2,69 Mg m' 3 , para 0,2-0, 4 
m de profundidade. Provavelmente, o manejo da area, onde nao se incorpora grande 
quantidade de materia organica, pode estar contribuindo para essa uniformidade. 

A densidade do solo (Figura 23a e 24a) apresentou uma relagao inversamente 
proporcional para as variaveis porosidade total (Figura 23b e 24b) e macroporosidade 
(figura 23c e 24c) e diretamente proporcional para a variavel microporosidade (Figura 
25a e 26a). Evidentemente, com o aumento da densidade, suas partfculas tendem a 
aproximar-se umas das outras, diminuindo assim o espago poroso do solo. No espa?o 
poroso que mais e afetado por esta aproximagao sao os poros de aeragao, e isso foi 
observado neste estudo para relagao entre densidade do solo e porosidade total, para 
ambas profundidades. 

O mesmo que ocorreu entre a densidade e a porosidade total do solo pode ser 
aplicado para a macroporosidade, ja que ha tambem uma correlagao negativa entre 
ambas. Percebe-se na Figura 23a e 23c, tambem na Figura 24a e 24c que os picos da 
densidade correspondem a depressao da macroporosidade. O aumento da densidade 
produz uma diminuigao inicial dos poros maiores. Na area de estudo, isto esta associado 
ao uso de implementos agrfcolas, principalmente arado e grade e ao trafego de 
maquinarias agrfcolas na area. 

O mapa de superffcie da Figura 23a mostra ainda uma grande depressao nas 
coordenadas 18 m x 20 m para densidade do solo, na profundidade de 0-0,2 m. 
Comportamento inverso foi observado para porosidade total e macroporosidade 
(Figuras 23b e 23c); Isto se deu devido a ocorrencia de um formigueiro exatamente 
nesta parte do terreno. O mapa da microporosidade se comportou de maneira identica ao 
da densidade do solo, ja que ambos tem uma relagao positiva (Figura 25a). Esta relagao 
pode estar associada ao manejo do solo e aos implementos que nele sao utilizados. 
Outro fator, citado anteriormente, e o fato da area nao ter nenhum tipo de cobertura 
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vegetal, ficando desprotegido contra os impactos das gotas de chuva, vento, dentre 
outros fatores. 

Na Figura 23a, observa-se que a cor rubra sobressai a cor negra, em ambas 
profundidades, sendo um possfvel indicador de compactagao do solo, pois o mapa 
aponta valores relativamente altos, mostrando a necessidade de seu controle com um 
manejo mais adequado. 

Quando se avalia os mapas de superffcie relacionados com a umidade do solo 
obtida a uma tensao de 0,033 MPa (Figura 25b e Figura 26b) e 1,5 MPa (Figura 25c e 
26c), percebe-se uma relagao positiva com a microporosidade. Para camada superficial 
e subsuperficial, os maiores valores, para as referidas variaveis se encontraram na 
porgao Norte da area. 

A textura do solo possui uma relagao com a disponibilidade de agua para as 
plantas. Assim, o teor de argila e um fator determinante para esta disponibilidade, pela 
maior area superficial de suas partfculas, podem reter maior volume de agua. Como 
citado anteriormente, o argilomineral predominante seja a caulinita, que tem cargas 
neutras, porem pode ocorrer cargas negativas na superffcie basal. Isso explica a 
interagao entre o teor de argila (Figura 19c e 20c) e a umidade obtida nas tensoes 0,033 
MPa e 1,5 MPa (Figura 25b, Figura 26b e Figura 25c, 26c). 

Verificou-se que a umidade obtida nas tensoes 0,033 e 1,5 MPa se comportaram 
de forma parecida com o teor de argila para as duas profundidades. Apresentando os 
maiores valores no final da area, nas coordenadas 20 m x 12 m, para profundidade de 
0,2-0, 4 m. 

Nota-se que os mapas das variaveis agregados 4,76-2,00 mm, 2,00-1,00 mm, 1,00- 
0,50 mm, 0,50-0,25 mm e <0,25 mm, apresentaram valores entre 100 e 900, 0 e 130, 0 e 
130, 0 e 120 e 100 e 650 g kg agregados' 1 , respectivamente (Figuras 27abc e 29ab), para 
profundidade de 0-0,2 m. Na profundidade de 0,2-0, 4 m, os valores variaram entre 250 
e 850 g kg agregados' 1 , 10 e 140 g kg agregados' 1 , 10 e 180 g kg agregados' 1 , 10 e 130 g 
kg agregados' 1 e 100 e 440 g kg agregados' 1 , respectivamente. 

Os agregados com diametro 4,76-2,00 mm predominaram em ambas 
profundidades. No entanto, houve um aumento com o incremento de profundidade, o 
que se constatou para os agregados com diametro 2,00-1,00 mm, 1,00-0,5 mm e 0,5- 
0,25 mm. Aqueles agregados com diametro menor que 0,25 mm diminuiu com o 
aumento da profundidade (Figura 27, 28, 29 e 30). 
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Figura 19 Distribui^ao espacial das variaveis areia, silte e argila total na 
profundidade de 0-0,2 m. 
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Figura 20. Distribui^ao espacial das variaveis areia, silte e argila total na 
profundidade de 0,2-0,4 m. 
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(a) 




Figura 21. Distribui^ao espacial das variaveis argila dispersa em agua, grau de 
floculacao e densidade das particulas na profundidade de 0-0,2 m. 
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Figura 22. Distribui^ao espacial das variaveis argila dispersa em agua, grau de 
flocula^ao e densidade das particulas na profundidade de 0,2-0, 4 m. 
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Figura 23. Distribui^ao espacial das variaveis densidade do solo, porosidade 
total e macroporosidade na profundidade de 0-0,2 m. 
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Figura 24. Distribui^ao espacial das variaveis densidade do solo, porosidade 
total e macroporosidade na profundidade de 0,2-0, 4 m. 
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Figura 25. Distribui^ao espacial das variaveis microporosidade, umidade 
0,033 MPa e umidade 1,5 MPa na profundidade de 0-0,2 m. 
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Figura 26. Distribui^ao espacial das variaveis microporosidade, umidade 
0,033 MPa e umidade 1,5 MPa na profundidade de 0,2-0,4m. 
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Figura 27. Distribui^ao espacial das variaveis estabilidade de agregados 4,76-2,00 
mm, estabilidade de agregados 2,00-1,00 mm e estabilidade de agregados 
1,00-0,50 mm na profundidade de 0-0,2 m. 
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Figura 28. Distribui^ao espacial das variaveis estabilidade de agregados 4,76-2,00 
mm, estabilidade de agregados 2,00-1,00 mm e estabilidade de agregados 
1,00-0,50 mm na profundidade de 0,2-0, 4 m. 
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Figura 29. Distribui^ao espacial das variaveis estabilidade de agregados 0,50-0,25 
mm e estabilidade de agregados menores que 0,25 mm na profundidade de 
0-0,2 m. 
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Figura 30. Distribui^ao espacial das variaveis estabilidade de agregados 0,50-0,25 
mm e estabilidade de agregados menores que 0,25 mm na profundidade de 
0,2-0, 4 m. 
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5. CONCLUSOES 

Com base nos resultados obtidos neste estudo, pode-se concluir que: 

1. As variaveis grau de floculagao, densidade do solo, na camada de 0-0,2 m e 
densidade das partfculas na camada de 0,2-0, 4 m apresentaram um efeito pepita 
relativo de 100% (efeito pepita puro), havendo independencia entre os pontos 
amostrados. 

2. Houve predommio do grau de dependencia espacial forte para as variaveis 
analisadas, tanto para aqueles influenciados pelos fatores de formagao do solo, 
como aqueles influenciados por fatores externos. 

3. O manejo aplicado ao solo diminuiu a dependencia espacial das caracterfsticas 
intrfnsecas como areia, silte e argila, consideradas pouco influenciadas pelo 
manejo. Por outro lado, propriedades altamente influenciadas pelo manejo como 
estabilidade de agregados e umidade na capacidade de campo apresentaram uma 
dependencia espacial significativamente aumentada. 

4. O manejo interferiu de modos diferenciados nos atributos analisados. Conferiu 
maior homogeneidade nas propriedades estruturais, implicando em menor 
numero de amostras para representar a area. Entretanto, em relagao as 
caracterfsticas texturais, o manejo aumentou a anisotropia implicando na 
necessidade de coletar maior numero de amostras para representar a area. 
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5. Os dados do alcance, para as variaveis analisadas, devem ser utilizados como 
base para delineamentos experimentais mais adequados e tambem como 
importante informagao acerca do manejo especffico local. 
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